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1. Introduccion 

La logica borrosa, o difusa, se puede considerar una disciplina dentro de la Inteligencia 
Artificial que permite representar el conocimiento acerca de un dominio y realizar procesos 
de inferencia o razonamiento. Para ello, proporciona un fonnalismo o lenguaje formal de 
representacion del conocimiento basado en la logica y las matematicas. 

La logica borrosa permite representar matematicamente conceptos que pertenecen a 
conjuntos cuyas fronteras no se encuentran nltidamente definidas, sino que estas son vagas, 
borrosas o difusas. Asl, por ejemplo, se puede establecer, con total seguridad, que una 
persona que mide 220 cm. pertenece al conjunto de las personas altas. Pero, £donde esta la 
frontera entre las personas consideradas altas y las personas de estatura media?. Se podrla 
establecer en este caso como valor frontera h = 180 cm.; de tal fonna que las personas cuya 
estatura sea igual o superior a 180 pertenecen al conjunto de las personas altas. Entonces, 
£una persona que mida 179 cm. no es alta?. No parece que se ajuste a la realidad. De hecho, 
no existe esa frontera. Cualquier punto que se elija resultara inadecuado en los valores de su 
entomo. En contraposicion con la teorla clasica de conjuntos, donde los elementos del 
dominio pertenecen o no pertenecen a un conjunto en base a unas reglas que establecen una 
frontera nitida, la pertenencia o no a un conjunto difuso o borroso resulta vaga o imprecisa 
en la frontera. 

La figura 1.1 muestra la diferencia existente entre frontera nltida e imprecisa. La figura 
1.1a) muestra la recta real que representa valores de temperaturas corporales (T), medidas en 
grados centlgrados (°C). El conjunto de las personas que tienen fiebre se establece mediante 
la expresion T > 37 °C. Por tanto, este dominio asl representado establece una frontera nltida 
(representada mediante una barra vertical) entre las personas que tienen fiebre y las que no, lo 
cual, se puede considerar, se ajusta bastante bien a la realidad. La figura 1.1b) muestra otra 
recta real que representa valores de altura de personas (h), medida en centimetros. En este 
caso no existe una frontera nitida (barra vertical) que separe a las personas altas del resto. Por 
el contrario, existe un tramo de indefinicion (vaguedad), en donde la pertenencia o no al 
conjunto de las personas altas no esta clara. Este tramo se establece en el intervalo h E 
(170,180). Por tanto, las personas que miden menos de 170 cm. no son altas, las personas que 
miden mas de 180 cm. si son consideradas asi, mientras que existe una frontera imprecisa de 
pertenencia en el intervalo (170,180) cm. Efectivamente, en este caso, es razonable pensar 
que esta forma de representar el conocimiento se ajusta mas a la realidad, a la forma que 
tienen las personas de percibir la realidad. Se da la paradoja, en este caso, de que la 
imprecision es mas precisa. Suponiendo que la frontera nitida establecida para la fiebre en 
la figura 1.1 a) es correcta, se extrae tambien la conclusion de que la logica borrosa no 
representa ni mejor ni peor el conocimiento que se tiene acerca de un dominio, sino que es un 
lenguaje fonnal de representacion adicional que puede ser utilizado para la construccion de 
sistemas inteligentes mas ajustados a la realidad. Debe ser una eleccion del ingeniero del 
conocimiento, infonnatico o disenador del sistema inteligente elegir el tipo de representacion 
mas adecuado al dominio o problema a resolver. 
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Figura 1.1. Fronteras nitidas e imprecisas. 


a. Razonamiento aproximado 

El razonamiento aproximado es una capacidad del ser humano por la cual es capaz de 
obtener conclusiones utiles a partir de conocimiento incomplete, impreciso o con cierto grado 
de incertidumbre. La logica tradicional se basa en que las clausulas, reglas o afirmaciones son 
ciertas o falsas con caracter absolute. Por ello, si bien esta es una forma valida de 
representacion del conocimiento y razonamiento en determinados dominios, la inteligencia 
artificial emplea modelos matematicos adicionales capaces de simular la capacidad humana 
para realizar razonamiento aproximado. Entre ellos, se encuentran la logica difusa, que 
proporciona un mecanismo matematico para realizar procesos de razonamiento con 
imprecision; y las redes bayesianas, o los sistemas MYCIN y Prospector, que permiten 
simular procesos de razonamiento con incertidumbre. 

Se considera que el conocimiento es impreciso, borroso o difuso, cuando las fronteras 
correspondientes a una determinada categoria son vagas. Esto es, hay un cierto numero de 
elementos para los que no se puede asegurar su pertenencia mediante valores dicotomicos: 
verdadero y falso o, cero y uno; sino que hay grados de pertenencia parcial. Asi por ejemplo, 
para la afirmacion “el dia es caluroso”, existen ciertos valores de temperatura ambiente para 
los que no se puede afirmar o rechazar dicha afirmacion. Diferente es el conocimiento con 
incertidumbre. La regia “Si fiebre, entonces gripe” es falsa si no se anade a este conocimiento 
cierto grado de incertidumbre, por ejemplo, mediante una medida de probabilidad (p). De 
esta forma, la regia anterior queda como: “Si fiebre, entonces gripe; p = 0.1”, siendo, en este 
caso, verdadera 1 . 

Con el fin de aclarar la diferencia entre imprecision e incertidumbre, se plantea el 
siguiente ejemplo: suponga el suceso aleatorio de lanzar un dado al aire, ( ;,quc numero 
saldra?. Cualquiera de los seis posibles valores es una afirmacion falsa: “saldra el 6”. Por 
tanto, es necesario emplear una forma de representar el conocimiento de forma aproximada. 
Utilizar en este caso la imprecision, implicaria que se podrian dar situaciones en las que no 
queda perfectamente determinado que numero sale una vez que el dado ha sido lanzado. Se 
estaria representando, por ejemplo, la posibilidad de que el dado quedara sobre una de sus 
aristas. Dado que este no es el caso y que lo que se desea es representar el conocimiento 
acerca del futuro (que numero saldra), es necesario emplear una forma de representar el 
conocimiento con incertidumbre: “saldra el 6; p = 1/6”, lo cual es estadisticamente cierto. Por 
tanto, se tiene que cada uno de los seis sucesos de variable aleatoria son precisos, pero existe 


1 El valor de probabilidad es inventado, para ilustrar el ejemplo. 
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incertidumbre acerca del numero que saldra. De hecho, en este caso, el conocimiento que 
realmente se esta representando es la ignorancia que se tiene hacia un hecho. Visto asl, la 
ignorancia es conocimiento, o mejor, una fonna de simular inteligencia. Efectivamente, tener 
la capacidad de responder “no lo se” es una fonna de exhibir inteligencia, dado que es mejor 
que dar una respuesta enonea. El mecanismo que tiene el uso de medidas de probabilidad 
para representar la ignorancia que se tiene acerca de un hecho es la asignacion de los mismos 
valores de probabilidad a todos los sucesos de la variable aleatoria: equiprobabilidad en todos 
los sucesos. Asl, en el caso del ejemplo del dado, la contestacion a la pregunta: “«^que numero 
saldra?”; la respuesta tambien puede ser: “no lo se”; lo que equivale a asignar un valor de 
probabilidad de 1/6 a cada uno de los 6 posibles resultados que se pueden obtener. 

Las redes bayesianas emplean, con gran exito, medidas de probabilidad condicionada 
para representar el conocimiento con incertidumbre. Un inconveniente que puede suscitar 
esta tecnica es la dilicultad de diferenciar entre ignorancia total acerca de un hecho, e igual 
certeza para todos los sucesos de una variable aleatoria. Suponga que se disena un sistema 
inteligente para el diagnostico de enfennedades vlricas. Una de tales enfermedades podrla ser 
la gripe, por lo que se declara la variable “gripe” con dos posibles valores (nltidos): “si” y 
“no”. Por tanto: gripe E {si, no}. Suponga que para diagnosticar la enfermedad, se lanza una 
moneda al aire. Si se obtiene cara, la respuesta es que el paciente tiene gripe. Si sale cruz, no 
la tiene. En este caso, la probabilidad de tener la enfennedad es 0.5, por tanto, gripe=“si”; 
p=0.5 es cierto. Tambien ocurre que gripe=“no”; p=0.5 es cierto. Evidentemente, dado el 
sistema de “diagnostico” empleado, se esta representando que se tiene ignorancia acerca del 
hecho de padecer la enfermedad. Suponga ahora la situacion en la que un paciente entra en 
una consulta medica porque tiene sospechas de gripe. El medico entonces hace una 
exploracion y encuentra signos de la enfennedad. Concluye, tras un proceso de inferencia en 
base a los hallazgos exploratorios, que en el caso concreto de ese paciente, la probabilidad de 
padecer gripe es del 50%. Se tiene, en este caso, con el mismo valor de probabilidad que 
anteriormente, certeza a favor del hecho de padecer la enfennedad. 

Para evitar este problema de diferenciacion entre ignorancia e igual certeza ante 
diferentes sucesos, el sistema inteligente MYCIN [Buchanan, 1984] emplea el concepto de 
factor de certeza (CF, en ingles). El factor de certeza que se asigna a un hecho o hipotesis 
toma valores entre -1 y 1: CF=[-1,1]. Asl, dada una hipotesis h, CF (h) = -1 indica que se 
esta totalmente seguro que no h; CF (h) = +1 indica que se esta seguro que se cumple h; los 
valores intermedios modelan diferentes grados de certeza a favor (positivos) o en contra 
(negativos) del hecho o hipotesis. Se reserva, en este caso, el valor cero, CF (h) = 0, para 
indicar ignorancia acerca de un hecho. 

Por tanto, y como conclusion, la representacion del conocimiento con imprecision 
pennite modelar dominios en los que se presentan categorlas o conjuntos cuyas fronteras no 
estan claramente definidas. Por su parte, la representacion del conocimiento con 
incertidumbre pennite modelar dominios en donde las afirmaciones no son estrictamente 
ciertas o falsas en todos los casos, aunque los valores que toman las variables son precisos. 
La incertidumbre pennite tambien representar la ignorancia total acerca de un hecho. 
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2.Fundamentos de la logica borrosa 

La logica borrosa proporciona un modelo matematico-logico para la construccion de 
sistemas inteligentes que manejan conocimiento con imprecision. El termino sistema 
inteligente es muy general, y se emplea para hacer referenda a cualquier sistema que 
desarrolla procesos que se pueden encuadrar dentro del area de la inteligencia artificial. 
Dentro de los sistemas inteligentes, se encuentran los sistemas basados en el conocimiento, 
cuya caracteristica principal es la representacion explicita del conocimiento mediante un 
lenguaje formal de representacion, que se emplea para llevar a cabo una tarea considerada 
como inteligente : diagnostico, planificacion o toma de decisiones entre otras. La figura 2.1 
muestra el esquema general de un sistema basado en el conocimiento, formado por una base 
de conocimiento y un motor de inferencia. Un sistema basado en el conocimiento recibe 
como entrada un conjunto de hechos que se cumplen para un caso particular y la pregunta 
que se formula al sistema en base a los hechos. Se dispara entonces un proceso de inferencia 
mediante su motor de inferencia, utilizando para ello la base de conocimiento, y se devuelve 
un resultado o respuesta. 


Hechos 



Pregunta que se fonnula 


Base de conocimiento 


Motor de inferencia 


Respuesta 



Figura 2.1. Estructura general de un sistema basado en el conocimiento. 

La base de conocimiento emplea un formalismo de representacion declarativa (no 
algoritmico) de los conceptos de un dominio especifico, sus propiedades, relaciones y 
elementos individuales que puedan aparecer. Existen diferentes formalismos de 
representacion del conocimiento, los cuales, en base a sus caracteristicas particulares, pueden 
resultar mas adecuados para representar los conocimientos que se manejan en uno o otro 
dominio especifico, como por ejemplo, la incertidumbre o la nitidez. La logica borrosa, por 
su parte, proporciona un lenguaje formal de representacion del conocimiento con 
imprecision. 

La logica borrosa pennite la definicion de un conjunto de variables que maneja el 
sistema inteligente: significado, unidades de medida, dominio de valores numerico y dominio 
de valores cualitativos o etiquetas linguisticas. Es precisamente en la definicion de este 
ultimo dominio de valores donde se encuentra la diferencia mas significativa de los sistemas 
inteligentes basados en logica borrosa. Suponga que se desea manejar y, por tanto, se define 
la variable T, temperatura ambiente en Madrid, medida en grados centigrados (°C) y cuyo 
dominio de valores es T E [-10,50] °C. Gracias al uso de la logica borrosa, podemos emplear 
tenninos o categorias para definir temperaturas cuyas fronteras no estan bien definidas, como 
por ejemplo: frio, templado, caluroso o alrededor de 40 °C (—40). Por tanto, ademas del 
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dominio de valores numerico definido para T, se define tambien el conjunto de etiquetas 
lingiiisticas o valores cualitativos, cuyos valores numericos no estan claramente dcflnidos: 
T={frlo, templado, caluroso, -40}. 

Tomando, por ejemplo, el valor T = caluroso, es facil darse cuenta que representa una 
categorla o un conjunto de valores numericos de temperaturas cuya frontera no se puede 
establecer en ningun valor concreto. Se podrla asegurar que los valores T E [30,50] 
pertenecen al conjunto caluroso. Si se acepta que T = 30 °C es una temperatura que se puede 
considerar como calurosa, entonces no es razonable suponer que T = 29 °C no lo sea. Por 
tanto, existe una frontera difusa, borrosa o vaga entre los valores de temperatura que se 
consideran calurosos y los que no. La figura 2.2 representa esta situacion. 



frontera 

borrosa 

caluroso 





-10 0 +10 +20 +30 +40 +50 T °C 


Figura 2.2. El conjunto de valores de temperatura (T) que pertenecen a la etiqueta 

T=caluroso tiene una frontera difusa. 

La base de conocimiento de un sistema inteligente basado en logica borrosa incluye un 
conjunto de reglas de la forma: 

van = vali OPRi van = vah OPR 2 ... OPR n -i var n = val n , —» varo = valo 

donde varo, van, ..., var n es un conjunto de variables que se han definido anteriormente; valo, 
vali, ..., vain es el conjunto de valores (numericos 0 cualitativos) que pueden tomar sus 
correspondientes variables; OPRi, OPR 2 , ..., OPR n -i son los operadores logicos de 
disyuncion (OR) 0 conjuncion (AND), y —» representa el operador de implicacion. Asi mismo 
es posible utilizar el operador logico unario de negacion de la forma: ^ (van = vali). 

El motor de inferencia de un sistema basado en el conocimiento es especifico para el 
formalismo de representacion del conocimiento empleado, y es independiente del dominio de 
aplicacion. Debe ser capaz de realizar procesos de inferencia con cualquier conjunto de 
conocimientos representados por el fonnalismo elegido. En el caso de la logica borrosa, el 
mas comunmente empleado es la regia composicional de inferencia 0 modus ponens 
generalizado, siendo el motor de inferencia de un controlador borroso un caso particular de 
este. 


a. Gonjuntos borrosos 

Se define un conjunto borroso como aquel en donde la pertenencia de sus elementos se 
define mediante la siguiente funcion de pertenencia, funcion de distribucion de 
posibilidad 0 de forma mas simple, funcion de posibilidad: 
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Siendo A el conjunto borroso, etiqueta linguistica o valor cualitativo perteneciente al 
dominio de valores cualitativos de la variable x; y .T cs el dominio de valores numericos de 
esta misma variable. La funcion de posibilidad pa (x) asigna a cada valor de x G SC un valor 

de posibilidad (y no probabilidad) de pertenencia al conjunto borroso A, comprendido entre 
cero (no pertenece) y uno (si pertenece). Los valores de posibilidad intermedios representan 
una pertenencia parcial de los valores numericos de x al conjunto A. Es importante subrayar 
que la funcion de pertenencia asigna valores de posibilidad de pertenencia de un valor 
numerico a un conjunto. Posibilidad se puede definir como la aptitud, potencia u ocasion para 
que algo exista o suceda (WordReference, 2015), en este caso, la pertenencia a un conjunto. 
No confundir este tennino con probabilidad: Indice por el cual se mide la frecuencia de un 
acontecimiento determinado mediante la realizacion de un experimento aleatorio del que se 
conocen todos los resultados (Wikipedia, 2015). Poniendo un ejemplo que aclare esta 
distincion, se podrla afirmar que “probablemente llueva en invierno en Madrid”, incluso con 
una cierta medida de probabilidad, porque se tiene experiencia pasada de que as! ha sucedido: 
existe una cierta frecuencia del acontecimiento. Por el contrario, serla incorrecto afirmar que 
“probablemente conocere al Presidente de los Estados Unidos”, aunque si es posible porque 
existe la potencia para ello, habida cuenta de que ambos estamos vivos. 

Para el caso del ejemplo mostrado en la figura 2.2, se podria definir el conjunto borroso 
“caluroso” mediante la siguiente funcion de distribucion de posibilidad, representada 
graficamente en la figura 2.3: 

Pcaiuroso: T a -> [0,1]; T a = [-10,50] °C. 

T * gcaluroso (T) 



Figura 2.3. Funcion de distribucion de posibilidad de pcaiuroso (T). 

En este caso se tiene por seguro que T = 40 °C es caluroso, y asi pcaiuroso (40) = 1. 
Tambien se sabe que T = 0 °C no pertenece al conjunto caluroso, por lo que pcaiuroso (0) = 0. 
Para los valores de T cuya pertenencia al conjunto “caluroso” es imprecisa, la funcion 
pcaiuroso(T) devuelve valores intermedios entre cero y uno. Por ejemplo pcaiuroso (23) = 0.3. La 
definicion del conjunto borroso “caluroso” mediante la funcion de pertenencia pcaiuroso (T) 
permite responder adecuadamente a la pregunta anteriormente planteada: si se considera que 
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una temperatura de 30 °C es calurosa, £se puede suponer que T = 29 °C tambien lo es?. En 
este caso, la respuesta es si, aunque con una posibilidad inferior, 0.9, porque p C aiuroso(29)= 0.9. 

Las funciones de posibilidad pueden tener, en general, cualquier forma, aunque se suele 
op tar por una de las tres mostradas en la figura 2.4: trapezoidal, en el caso de ptempiado(T); 
triangular, en el caso de p~ 2 o (T); o nitida como ocurre con p4o (T). Se puede considerar que 
Pcaiuroso (T), mostrada en la figura 2.3, tiene una forma trapezoidal. No se representa el lado 
derecho por estar fuera del dominio de valores definidos para las temperaturas en Madrid, 
que, se recuerda, es [-10,50] °C. 



Figura 2.4. Tres formas diferentes para las distribuciones de posibilidad: trapezoidal, 

triangular y nitida. 

La distribucion de posibilidad p ~20 (T) (leido como mu sub alrededor de veinte ) define 
el conjunto borroso formado por todos los valores de temperatura que se encuentran 
alrededor de T = 20. En general, para representar conjuntos borrosos formados por valores 
que se encuentran en el entomo de uno central, se emplean distribuciones de posibilidad cuyo 
centro de gravedad se situa en dicho valor central. En el caso de p~ 2 o (T) se ha optado por 
una distribucion de posibilidad triangular, con centro de gravedad en T = 20, aunque tambien 
hubiera sido valida una distribucion de posibilidad trapezoidal. En este caso, se ha 
representado una cantidad aproximada, pero se dan situaciones similares en los casos de 
querer representar que alguien o algo es puntual (en funcion del tiempo medido en minutos, 
t), que una trayectoria recta (en funcion del angulo en radianes, ca) o que un sistema 
(controlador) comete un error nulo (en funcion del error, e). En todos ellos, se podria 
representar la puntualidad, la trayectoria o el error en una recta real. Los valores a la 
izquierda del centro de gravedad indicarian, respectivamente, antes de tiempo, trayectoria 
desviada hacia la izquierda y error negativo, por ejemplo. Los valores a la derecha del centro 
de gravedad representarian las situaciones contrarias: tarde, desvio a la derecha y error 
positivo. Para representar las situaciones en_punto, trayectoria recta y error nulo, se 
emplearian funciones de distribucion de posibilidad triangulares o trapezoidales cuyo centro 
de gravedad estuviera situado en el valor nitido representado por la etiqueta linguistica: 0, n/2 

y 0, respectivamente. La figura 2.5 muestra ejemplos de distribuciones de posibilidad para 
estos tres ejemplos. 
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Figura 2.5. Funciones de distribucion de posibilidad para representar valores 
pertenecientes al entorno de uno central: alguien o algo es puntual, p e n_ P unto (t), una 
trayectoria recta, p rec to ( 10 ), y un sistema de control que comete un error nulo, p nu io (e). 


Volviendo a la figura 2.4, la funcion de distribucion de posibilidad p4o (T), representa el 
valor nitido 40. Efectivamente, los conjuntos borrosos permiten representar tambien 
conjuntos o valores nitidos, que tienen fronteras precisas. En el caso de p 4 o (T), cualquier 
valor en el entorno de 40 tiene posibilidad cero de pertenecer a T = 40: p 40 (40 ± s) = 0, s—>0; 
mientras que p4o (40) = 1. Matematicamente, p 4 o(T) se formula como: 



1, T= 40 
0, en otro caso 


Se muestra a continuacion un ejemplo, dentro del dominio de las temperaturas 
ambientales, de definicion de variables, dominios de valores numericos y cualitativos y 
funciones de distribucion de posibilidad de un sistema inteligente basado en logica borrosa 
que maneja unicamente la variable T. 


1. Definicion de las variables del sistema inteligente, significado y unidades de medida. 
• T: temperatura ambiente en Madrid, medida en grados centigrados (°C). 


2. Definicion de los dominios de valores numericos y cualitativos de las variables. 

. TG [-10,50] °C. 

• T = {frio, templado, caluroso, alrededor de 20, 40} 

3. Definicion de las funciones de distribucion de posibilidad asociadas a los valores 
cualitativos de cada una de las variables definidas (en este caso T). Ver figura 2.6. 



Figura 2.6. Representacion grafica de las funciones de distribucion de posibilidad 
asociadas a los valores cualitativos definidos para la variable T. 
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El elemento caracterlstico de esta definicion de variables en un sistema inteligente basado 
en logica borrosa es el diseiio de las funciones de distribucion de posibilidad asociadas a las 
etiquetas linguisticas, conjuntos borrosos o valores cualitativos de a cada una de las variables. 

Es posible tambien definir conjuntos borrosos en n dimensiones mediante funciones 
de distribucion de posibilidad de la forma p p (xi, x n ). La figura 2.7 muestra un ejemplo 
que representa una funcion de distribucion de posibilidad para el conjunto borroso de las 
personas corpulentas en funcion de la altura (h), medida en cm., y el peso (p), medido en Kg, 

notado COmO pcorpulento (h, p). 



Figura 2.7. Ejemplo de funcion de distribucion de posibilidad en dos dimensiones: 

Pcorpulento (h, p). 


b. Operaciones con conjuntos borrosos 

Dadas dos distribuciones de posibilidad p P (x) y p q (x) correspondientes a los conjuntos 
borrosos p y q, es posible realizar las operaciones de interseccion (D) y union de conjuntos 
borrosos (U), asi como el calculo del conjunto complementario a uno dado (-). Estas 
operaciones se realizan a traves de la funcion de distribucion de posibilidad asociada a cada 
conjunto borroso (composicion de funciones). Suponga que las distribuciones de posibilidad 
p P (x) y p q (x) se corresponden con las de la figura 2.8. 



Figura 2.8. Representacion grafica de las distribuciones de posibilidad p p (x) y p q (x) 

A partir de la teoria clasica de conjuntos, es posible obtener, de forma intuitiva, el 
resultado de aplicar una operacion de interseccion y union entre los conjuntos borrosos p y q 
a traves de sus correspondientes funciones de distribucion de posibilidad, p P (x) y p q (x); asi 
como el complementario p, lo cual se muestra en la figura 2.9. 
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Figura 2.9. Resultado de aplicar las operaciones de intersection y union de p y q a traves 
de sus funciones de distribucion de posibilidad p p (x) y p q (x), y el complementario de p. 

Las operaciones con conjuntos borrosos (interseccion, union y complementario) tienen 
sus correspondientes operadores logicos: conjuncion (a), disyuncion (v) y negation 
respectivamente. A partir de la cuantificacion de la posibilidad de que sean ciertas “x es p” y 
“x es q” mediante p p (x) y p q (x), respectivamente, es posible expresar que “x es p y q” como 
p P Aq (x), “x es p 6 q” como p P v q (x) y “x es no p” como p- p (x). Asi mismo, se define el 
operador logico implicacion (—>■), que permite expresar que “si x es p, entonces x es q” como 
p P —> q (x). En este caso, resulta mas conveniente nombrar a p y q como etiquetas linguisticas o 
valores cualitativos en lugar de como conjuntos borrosos. Dicho de otra manera, las etiquetas 
linguisticas o valores cualitativos p y q son representados como conjuntos borrosos mediante 
la definicion de sus correspondientes funciones de distribucion de posibilidad p P (x) y p q (x), 
respectivamente. 

De las misma manera, generalizando a n dimensiones, dadas las distribuciones de 
posibilidad p P (xi, X2, ..., x n ) y p q (xi, X2, ..., x n ), es posible expresar la conjuncion, 
disyuncion, negacion e implicacion de sus correspondientes etiquetas linguisticas (p y q) 
como p P Aq (xi, X2, ..., Xn), p p v q (xi, X2, ..., X n ), P^p (xi, X2, ..., X n ), |4p-> q (xi, X2, ..., X n ), 
respectivamente. 

Mediante el uso de los operadores logicos, es posible representar el conocimiento con 
imprecision mediante afirmaciones y reglas borrosas que relacionan etiquetas linguisticas. 
Asi, volviendo al ejemplo de las temperaturas ambientales, es posible expresar las 
afirmaciones “hoy el dia es frio o templado” como pfdovtempiado (T), “no hace calor” como 
p - 'caiuroso (T) o la regia “si la temperatura es de alrededor de 20 °C, entonces el dia es 
templado” como p~20-uem P iad<>(T). 


c. Extension cilindrica 

Solo es posible componer (operar) dos distribuciones de posibilidad utilizando los 
operadores de conjuncion, disyuncion o implicacion cuando ambas funciones se encuentran 
expresadas mediante la misma variable o conjunto de variables. En efecto, todos los ejemplos 
empleados hasta el momenta se encuentran expresados en funcion de la misma variable, x o 
T. Sin embargo, la capacidad de representacion del conocimiento que ello proporciona es 
muy limitada, dado que no podrian expresarse, por ejemplo, reglas del tipo “si la distancia 
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recorrida con el taxi es corta, entonces el precio es barato”. En este caso, se estan manejando 
dos variables, distancia y precio, relacionadas mediante la operacion implicacion. Mas 
formalmente, lo que se define es “si distancia = corta, entonces precio = barato”, 
representado como: pcorta (distancia)—>pbarato(precio). 

En general, dadas dos distribuciones de posibilidad p P (x) y p q (y), donde x e y pueden 
ser n-dimensionales, no es posible componerlas porque ambas deben estar definidas en 
funcion de las mismas variables. Para ello, dada una distribucion de posibilidad p P (xi, X2, ..., 
x n ), se define la operacion extension cilindrica de p p (xi, X2, ..., x n ) con y a la distribucion de 
posibilidad p p (xi, X2, ..., x n , y) / V y p P (xi, X2, ..., x n , y) = p p (xi, X2, ..., x n ). Es decir, se 
introduce la variable y en la distribucion de posibilidad original, con lo que aumenta en una 
dimension (se extiende cilmdricamente), pero los resultados de posibilidad que devuelve la 
funcion unicamente dependen del vector original (xi, X 2 , ..., x n ). 

De esta forma, para poder componer dos distribuciones generales de posibilidad p P (x) y 
p q (y) para calcular, por ejemplo, la conjuncion entre p y q, es necesario calcular, en primer 
lugar, la extension cilindrica de p p (x) con y para obtener p p (x, y). Posterionnente, se calcula 
la extension cilindrica de p q (y) con x para obtener p q (x, y), que es, evidentemente, 
equivalente a p q (y, x). Por ultimo, puesto que ya se tienen las dos distribuciones de 
posibilidad expresadas en funcion de las mismas variables (x,y), se aplica el operador 
correspondiente, en este caso la conjuncion, para obtener p pAq (x,y). Logicamente, el operador 
de negacion no tiene este problema, dado que es un operador unario, por lo que opera 
unicamente con una distribucion de posibilidad. 

Volviendo al ejemplo de las distribuciones de posibilidad p C orta (distancia) y pbarato 
(precio), para poder calcular la implicacion entre corta y barato, representado a traves de sus 
correspondientes funciones de distribucion de posibilidad como p C orta (distancia) —> 
pbarato(precio), es necesario calcular la extension cilindrica de p C orta (distancia) con precio para 
obtener p CO rta (distancia, precio). Asi mismo, se calcula la extension cilindrica de 
Pbarato (precio) con distancia para obtener pbarato (distancia, precio). Finalmente, se aplica el 
operador implicacion para obtener la distribucion de posibilidad resultado: 

pcorta ->baiato(distancia, precio). 

La figura 2.10 muestra las extensiones cilindricas de p CO rta (distancia) con precio, y 
Pbarato (precio) con distancia. Observe como, por ejemplo, en el caso de p C orta (distancia, 
precio), los valores de posibilidad devueltos unicamente dependen de la variable distancia. Es 
decir, p CO rta (distancia) = p CO rta (distancia, precio), V precio. La inclusion de la variable precio 
en la funcion tiene el efecto de aumentar el dominio en una dimension: de ser una funcion 
R—>[0,1], pasa a ser R 2 —>[0,1]. De la misma forma, observe en la figura 2.10 que pbarato 
(precio) = pbarato (distancia, precio), V distancia. 
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Figura 2.10. Representacion grafica de las extensiones cillndricas de p CO rta (distancia) con 

precio y pbarato (precio) con distancia. 


d. Funciones t-norma, t-conorma y negacion 

Una t-norma (norma triangular) es una funcion T definida como: 

T: [0, 1] x [0, 1] - [0, 1] 

(x, y) —► T (x, y) 

donde x representa el producto cartesiano. Una t-norma cumple, ademas, las siguientes 
propiedades, V x, y, z, t E[0, 1 ]: 

• Propiedad conmutativa: T (x, y) = T (y, x). 

• Propiedad asociativa: T (T (x, y), z) = T (x, (T (y, z)). 

• Elemento neutro, 1: T (x, 1) = x 

• Monotonia: x < z, y < t —>■ T (x,y) < T (z, t). 

Una t-norma representa la interseccion dentro de un reticulo. Se emplea en logica borrosa 
para calcular la interseccion de conjuntos borrosos o la conjuncion entre etiquetas 
linguisticas o valores cualitativos en el ambito de la logica. 

Sean las distribuciones generales de posibilidad p p (x) y p q (x), expresadas en funcion de 
la misma variable o vector x E3T, siendo SC el dominio de valores (o vectores) numericos de 

dimension n que puede tomar x. Si no ocurriera asi, se aplican las operaciones de extension 
cilindrica a ambas funciones para conseguir que queden expresadas en funcion de la misma x. 

p P (x) y p q (x) definen los conjuntos borrosos p y q mediante medidas de posibilidad de 
pertenencia de cada uno de los valores x E.T a p y q respectivamente: p P : SC —> [0, 1]; p q : SC 

—» [0, 1]. En el ambito de la logica, p p (x) y p q (x) representan las afirmaciones “x es p” y “x 
es q” respectivamente. Si se aplica la funcion t-norma a cada uno de los valores de 
posibilidad de p p (x) y p q (x) V x E SC, se obtiene otra distribucion de posibilidad, en funcion 
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de la misma variable x, que representa la interseccion de los conjuntos borrosos p H q, o la 
conjuncion en el ambito de la logica: “x es p y q”. Esta nueva funcion de distribucion de 
posibilidad se expresa formalmente como: 

p pAq (x): SC —> [0, 1], 

Sea M el conjunto de todas las funciones de distribucion de posibilidad expresadas en 
funcion de la variable o vector x. V p p (x), p q (x) G M, la funcion t-norma, T, aplicada a p p (x) 

y p q (x), V x ESC se expresa como: 

H-P ( X ), Pq (x) —> p P Aq (x) 

Asl, para calcular p p (x) a p q (x), se aplica la funcion t-norma: T (p p (x), p q (x)) = p pAq 
(x), V x 6 Las funciones t-norma aplicadas a distribuciones de posibilidad cumplen las 
mismas propiedades, V p p (x), p q (x), p r (x), p s (x) G M : 

• Propiedad conmutativa: T(p P (x), p q (x))=T(p q (x), p p (x)). 

• Propiedad asociativa: T[p p (x), T (p q (x), pr(x))]=T[T(p p (x), p q (x)), p r (x)]. 

. Elemento neutro: la funcion de posibilidad pi (x) = 1 V x ESC; T(p p (x), pi (x))=p p (x) 

• Monotonia: 

Si V x G SC, p P (x) < p r (x), p q (x) < p t (x) —> T(p p (x), p q (x)) < T(p r (x), p t (x)). 

Con el fin de simplificar la notacion, se expresan las funciones de distribucion de 
posibilidad p P (x), p q (x), p r (x), p s (x) G M como x, y, z, t respectivamente, de tal fonna que 
T (p P (x), p q (x)) queda como T (x, y). La funcion de distribucion de posibilidad pi (x) = 1 V 
x G3T, elemento neutro de la t-norma, se expresa como 1, y po (x) = 0 V x ESC, como 0. 

Ejemplos de funciones t-norma son las siguientes, V x G SC: 

• t-norma del minirno (o de Godel): T(x, y)= min (x, y). 

• t-norma del producto, P: T(x, y) = P(x, y) = x • y. 

• t-norma de Lukasiewicz, W: 

T(x, y) = W(x, y) = max (0, x + y - 1). 

• t-norma drastica, Z: 


T(x,y) = Z(x,y) 


x, si y = 1 

y, si x = 1 
0, en otro caso 
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Se puede demostrar de forma sencilla que los cuatro ejemplos de funciones t-norma 
anteriores cumplen las propiedades de toda t-norma. 

La funcion de distribucion de posibilidad que representa la interseccion (o conjuncion) 
entre p y q en la figura 2.9 muestra el resultado de aplicar la t-norma del mlnimo a las 
distribuciones de posibilidad definidas graficamente en la figura 2.8. 

Una t-conorma (o S-norma) es una funcion S (tambien notada como J.) deflnida como: 

S: [0, 1] x [0, 1] - [0, 1] 

(x, y) —* S (x, y) 

donde x representa el producto cartesiano. Una funcion t-conorma cumple, ademas, las 
siguientes propiedades, V x, y, z, t E[0, 1]: 

• Propiedad conmutativa: S (x, y) = S (y, x). 

• Propiedad asociativa: S (S (x, y), z) = S (x, (S (y, z)). 

• Elemento neutro, 0: S (x, 0) = x 

• Monotonia: x < z, y < t —S (x,y) < S (z, t). 

La t-cononna se emplea en logica borrosa para calcular la union de conjuntos borrosos o la 
disyuncion de etiquetas lingulsticas o valores cualitativos en el ambito de la logica. 

Sean las distribuciones generales de posibilidad p p (x) y p q (x), expresadas en funcion de 
la misma variable o vector x E3T. Si se aplica la funcion t-conorma a cada uno de los valores 

de posibilidad de p p (x) y p q (x), V x ElSC, se obtiene otra distribucion de posibilidad, en 

funcion de las mismas variables x, que representa la union de los conjuntos borrosos p U q, o 
la disyuncion en el ambito de la logica: “x es p 6 q”. Esta nueva funcion de distribucion de 
posibilidad se expresa formalmente como: 

E P v q (x): 3C -*■ [0, 1] 

V p p (x), p q (x) E M, la funcion t-conorma, S, aplicada a p p (x) y p q (x), V x E.T se 
expresa como: 

S: M 2 —>■ M 

M-P (x), Pq (x) —► Ppvq (x) 

Asi, para calcular p p (x) v p q (x), se aplica la funcion t-conorma: S (p p (x), p q (x)) = p pVq 
(x), V x E.T. De nuevo, para simplificar la notacion, se plantean las siguientes equivalencias: 

sean p p (x), p q (x), p r (x), p s (x) E M; p p (x) = x, p q (x) = y, p r (x) = z, p s (x) = t. De esta 
forma, T (p p (x), p q (x)) = T (x, y). Asi mismo, la funcion de distribucion de posibilidad 
pi (x) = 1 V x E3T se expresa como 1, y po (x) = 0 V x E3T, como 0. 

Las funciones t-nonna aplicadas a distribuciones de posibilidad cumplen las mismas 

propiedades V x, y, z, t E M : 
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• Propiedad conmutativa: S(x, y) = S(y, x). 

• Propiedad asociativa: S [x, S (y, z)] = S [S (x, y), z], 

• Element® neutro: la funcion de posibilidad po (x) = 0, V x G3T; S (x, 0) = x 

• Monotonia: x < z, y < t —* S (x, y) < S (z, t). 

Ejemplos de funciones t-conorma son las siguientes, V x 

• t-conorma del maximo (o de Godel): S (x, y)= max (x, y). 

• t-conorma de la suma-producto, P’: S(x, y) = P’(x, y) = x + y - x • y. 

• t-conorma de Lukasiewicz, W’: S(x, y) = W’(x, y) = min (1, x + y). 

• t-conorma drastica, Z’: 


S(-Ly) = Z'(x, /t q (x)) 


x, si y = 0 

y, si x = 0 

1, en otro caso 


Se puede demostrar de forma sencilla que los cuatro ejemplos de funciones t-conorma 
anteriores cumplen las propiedades de las t-conormas. 

La funcion de distribucion que representa la union (o disyuncion) entre p y q en la 
figura 2.9 muestra el resultado de aplicar la t-conorma del maximo a las distribuciones de 
posibilidad definidas graficamente en la figura 2.8. 

La funcion de negacion, N, se define como: 

N: [0,1] —>• [0,1] 

x —> N (x) 

La funcion de negacion N cumple las siguientes propiedades, V x, y G[0, 1]: 

• Condiciones frontera: N (0) = 1; N (1) = 0. 

• Inversion de monotonia: x < y —> N (x) > N (y). 

La negacion se emplea en logica borrosa para calcular el conjunto borroso complementario 
o la negacion en el ambito de la logica. 

Sea la distribucion general de posibilidad p P (x). Si se aplica la funcion negacion a cada 
uno de los valores de posibilidad de p P (x) V x ElSC, se obtiene otra distribucion de 

posibilidad que representa el conjunto borroso complementario de p, o la negacion en el 
ambito de la logica: “x es no p”. Esta nueva funcion de distribucion de posibilidad se expresa 
formalmente como: 

p-p(x): 3C —► [0, 1], 
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V p p (x) E M, la funcion negacion, N, aplicada a p p (x), V x El3C se expresa como: 

N \M^M 
p P (x) -»■ p p (x) 

La funcion de negacion, aplicada a funciones de distribucion de posibilidad, cumple las 
mismas propiedades. V x, y EM, representando po (x) = 0, V x E 3C como 0 y pi (x) = 1, V 

x E 3C como 1: 

• Condiciones frontera: N (0) = 1; N (1) = 0. 

• Inversion de monotonia: Si V x E<2T, x < y —>■ N (x) > N (y). 

Ejemplos de funciones de negacion son las siguientes, V x E 3C\ 

• N (x)= 1 - x. 

• N (x)= [1 - x 2 ] 1/2 . 

Se podria demostrar de forma sencilla que los dos ejemplos de funciones de negacion 
anteriores cumplen las propiedades. 

La distribucion de posibilidad que representa el complementario (o negacion) de p en la 
figura 2.9 muestra el resultado de aplicar la funcion de negacion N (x) = 1 - x a la 
distribucion de posibilidad p P (x), definida graficamente en la figura 2.8. 

Notese que, al contrario de lo que ocurre en la logica bivalente, la conjuncion entre una 
afirmacion y su negacion en logica borrosa, calculada mediante la t-norma de sus 
correspondientes funciones de distribucion de posibilidad, no tiene porque ser necesariamente 
falsa (valor cero de posibilidad) para todos los valores del dominio. Es decir: 

3 p p (x) E Ml T (p p (x), p-p (x)) ^ po (x); siendo po (x) = 0 V x E 3C . 



Figura 2.11. V x E(a,b) Car,T (p p (x), p-, p (x)) ^ 0. Calculo realizado con la t-norma del 

minimo. 

La figura 2.11 muestra un caso generico en el que la conjuncion entre las afirmaciones 
“x es p” y “x es no p”, calculado como T (p p (x), p- p (x)), con la t-norma del minimo, obtiene 
como resultado una distribucion de posibilidad en donde se cumple que V x E (a,b) C .T, T 
(p P (x), p^,(x))^0. 
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Paralelamente, en teoria de conjuntos borrosos, tambien se cumple que, dado un 
conjunto borroso p definido mediante la funcion de pertenencia p p (x), p fl “’p ^ 0. Siendo 0 
el conjunto vaclo. 

e. Dualidad 

Un sistema inteligente basado en logica borrosa puede utilizar cualquier t-norma, t- 
cononna o funcion de negacion, siempre y cuando se cumplan las condiciones de dualidad 
entre ellas. 

V x, y G M, se dice que una t-norma T (x, y) y una t-conorma S (x, y) son duales con 
respecto a la negacion N (x), si y solo si V x G 3T: 

• N (T (x, y)) = S ( N (x), N(y)). 

• N(S(x,y)) = T(N(x),N(y)). 

En relacion a las t-normas y t-conormas expuestas como ejemplo anteriormente en la 
seccion 2.d, son duales con respecto a la negacion N(x) = 1 - x, con 1 = pi (x) = 1, V x G3T, 
los siguientes pares: 

• < mm (x, y), max (x, y) >, Godel. 

• < P (x, y), P’ (x, y) >, Producto, Suma-Producto. 

• < W (x, y), W’ (x, y) >, Lukasiewicz. 

• < Z (x, y), Z’ (x, y) >, Drastica. 

Se demuestra a continuacion, a modo de ejemplo, que la t-norma del producto y la t- 
cononna de la suma-producto son duales con respecto a la negacion N(x) = 1 - x, V x G 3C. 

Las condiciones necesarias y suficientes para que se de la dualidad son las siguientes: 

1. N(P(x,y)) = P’(N(x),N(y)). 

2. N(P’(x,y)) = P(N(x),N(y)). 

Para la primera condicion, sustituyendo a N por su correspondiente expresion 1 - x, se tiene: 
l-(P(x,y)) = P’(l-x,l-y) 

Sustituyendo ahora P y P’ por sus correspondientes expresiones: 

1 - x-y = 1 - x + 1 - y -[(1 - x) - (1 - y)] 

Desarrollando la parte derecha de la igualdad: 

1 - x-y = l- x+ l- y-[l-x-y + x-y] 

Se elimina el signo menos anterior a los corchetes: 

1 - x-y = l- x+ l- y- l+ x + y- x-y 
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Se tiene x, -x; y, - y; 1 - 1; por tanto, se simplifica, quedando : 

1 - x-y = 1 - x-y 

Para la segunda condicion, se desarrolla de fonna similar: 

l-(x + y-xy) = (l-x)-(l-y) 

Desarrollando ambos lados de la igualdad, se tiene: 
l-x-y+x-y=l-x-y+x-y 

Dado que se cumplen las dos leyes de la dualidad y que son condiciones necesarias y 
suficientes para ello, queda demostrado que P (x,y) = x ■ y, y P’ (x,y) = x + y - x • y son 
duales con respecto a la negacion N (x) = 1 - x. 

Para construir un sistema inteligente que maneje conocimiento con imprecision 
utilizando la logica borrosa, es necesario que las funciones t-norma T y t-conorma S 
elegidas sean duales con respecto a una negacion N. A partir de estas tres funciones, que 
penniten calcular relaciones entre conjuntos borrosos o etiquetas linguisticas a traves de sus 
correspondientes distribuciones de posibilidad, es posible representar el conocimiento en base 
a afirmaciones que combinen conjunciones, disyunciones y negaciones. Por ejemplo: “x es p 
y no q” o “la temperatura ambiente es calurosa y no fria”. Estas afirmaciones penniten 
construir los antecedentes de las reglas borrosas que forman parte de la base de conocimiento 
de un sistema inteligente basado en reglas borrosas. 


f. Implicacion difusa 

Sean las distribuciones generales de posibilidad p p (x) = x y p q (x) = y, expresadas en 
funcion de la misma variable o vector x GST, siendo SC el conjunto de valores o vectores 

numericos de dimension n que puede tomar x. Si no ocurriera asi, se aplicarian las 
operaciones de extension cilindrica a ambas funciones para conseguir que queden expresadas 
en funcion de la misma x. 

Sea M el conjunto de todas las funciones de distribucion de posibilidad expresadas en 
funcion de la variable o vector x. V p p (x), p q (x) G M, la funcion implicacion, J, se define 
como: 

p P (x), p q (x) pL P > q (x) 

La funcion implicacion recibe como parametro dos distribuciones de posibilidad 
expresadas en funcion de la misma variable o vector, y devuelve, como resultado, otra 
distribucion en funcion de las mismas variables que representa la implicacion de las 
correspondientes etiquetas linguisticas o valores cualitativos en el ambito de la logica. Asi, 
para calcular p p (x) —> p q (x), se aplica la funcion implicacion: J (p p (x), p q (x)) = J (x, y) = 

Pp > q ( x )- 
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Dada una t-conorma S y una negation N, se puede obtener una expresion para la 
funcion implication J. A partir de la logica propositional, dadas dos proposiciones A y B, se 
sabe que A—>B = “AvB. Sustituyendo la negation por la funcion de negation borrosa N y la 
disyuncion por la t-conorma S, se tiene: J (x, y) = S (N (x), y). Tomando: N (x) = 1 - x, con po 
(x) = 0 = 0, pi (x) = 1 = 1, V x (E SC, es posible obtener una expresion para la implication 

difusa en base a la t-conorma S elegida: 

• Con S (x, y) = max (x, y), se tiene: 

J (x, y) = max (1-x, y). 

• Con S (x, y) = P’(x, y) = x + y - x • y, se tiene: 

J (x, y) = 1 - x + x • y. 

• Con S (x, y) = W’(x, y) = min (1, x + y), se tiene: 

J (x, y) = min (1, 1 - x + y). 

• Con: 


Se tiene: 


S(i,y) = Z'(i,yH 


x, siy = 0 

y, si x = 0 

1, en otro caso 


J(i,y) 


y, si x = 1 
1 - x, si y = 0 
1, en otro caso 


A modo de ejemplo ilustrativo, se indica a continuation como se deduce la funcion de 
implication J a partir de la t-conorma de la suma-producto, S (x, y) = P’(x, y) = x + y - x • y, 
y la negacionN(x) = 1 -x, VxG J: 

Se tiene que J (x, y) = S (N (x), y). 

Sustituyendo N por su valor: 

J (x,y) = S (1-x, y) 

Puesto que S (x, y) = P’ (x, y) = x + • y - x • y: 

J(x, y) = 1-x+ y-(1-x) -y 

Operando: 

J (x, y) = 1-x + y - y + x y 
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Simplificando: 


J (x, y) = 1 -x + xy 

Utilizando la notacion tradicional: 

J (up (x), p q (x)) = pi (x) - p P (x) + p p (x) • p q (x) 

Dadas las proposiciones A y B, a partir de las equivalencias en logica proposicional, 
tambien se tiene que AvB = Av(~AaB). Puesto que A—>B = ~AvB, es posible obtener otro 
conjunto de expresiones validas para la implicacion difusa con N (x) = 1 - x, con po (x) = 0 = 
0, pi (x) = 1 = 1, V x G SC, una t-norma T cualquiera y la correspondiente t-conorma S dual 
con respecto a la negacion N: 

• Con < mm (x, y), max (x, y) >, se tiene: 

J (x, y) = max (1 - x, min (x, y)). 

• Con < P (x, y), P’ (x, y) >, se tiene: 

J (x, y) = 1 - x + x 2 • y. 

• Con < W (x, y), W’ (x, y) >, se tiene: 

J (x, y) = max (1 - x, y). 

• Con < Z (x, y), Z’ (x, y) >, se tiene: 

y, si x = 1 
1 - i, si x, y J 1 1 
1 ,x # l,y = 1 

Por tanto, dado N (x) = 1 - x, con po(x) = 0 = 0, pi (x) = 1 = 1, V x G 3C, una t-norma T 
cualquiera y la correspondiente t-cononna S dual con respecto a la negacion N, es posible 
deducir y emplear dos funciones de implicacion J diferentes, siendo ambas correctas desde el 
punto de vista formal. 

Implementar un motor de inferencia para un sistema basado en el conocimiento con 
logica borrosa como formalismo de representacion del conocimiento implica elegir una t- 
nonna para la conjuncion, una t-conorma para la disyuncion, una negacion y una funcion de 
implicacion. La t-nonna y t-conorma elegidas deben ser duales con respecto a la negacion. 
As! mismo, en general (excepto en el caso de implementar un controlador borroso), se debe 
elegir una funcion de implicacion que se deduzca de las anteriores. 

En el caso de los controladores borrosos, se emplea siempre la funcion de implicacion 
de Mamdani [Mamdani, 1977]. Dadas dos funciones de distribucion de posibilidad p P (x), p q 
(x) G M, la funcion de implicacion de Mamdani se podria definir como, V x G 3T: 

J (p P (x), p q (x)) = min (p p (x), p q (x)). 


J(i,y) =\ 
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O, de forma equivalente: 


J (x, y) = min (x, y). 

En base a la definicion de las cuatro operaciones en logica borrosa (conjuncion, 
disyuncion, negacion e implicacion) es posible construir bases de reglas que manejan 
tenninos con valores cualitativos del tipo: “Si la cantidad de monoxido de carbono es elevada 
y la tendencia es positiva o la visibilidad no es clara, entonces la potencia de ventilacion debe 
ser alta”. Esta regia combina las cuatro operaciones logicas e involucra a cuatro variables 
diferentes (cantidad de monoxido de carbono, tendencia, visibilidad y potencia de 
ventilacion), cada una de las cuales con un valor cualitativo asignado (elevada, positiva, clara 
y alta, respectivamente). Tanto las variables involucradas como sus valores cualitativos y las 
funciones de distribucion de posibilidad correspondientes deben haber sido definidas 
previamente para la construccion de un supuesto sistema basado en el conocimiento que 
controle la ventilacion de un recinto. 


g. Ejemplos 

1. Sean las distribuciones generates de posibilidad ju p (x) y p q (y). iQue operaciones hay que 
realizarpara poder calcular la distribucion de ju p (x) —*■ p p (y)?. 

Para poder calcular la implicacion entre dos distribuciones de posibilidad es necesario que 
ambas esten definidas en funcion de las mismas variables. Para ello, en primer lugar, es 
necesario calcular las correspondientes extensiones cilindricas con x e y. Posterionnente, se 
calcula la implicacion utilizando la funcion J correspondiente. Por tanto: 

• Calculo de la extension cilindrica de p P (x) con y, para obtener p p (x, y). 

• Calculo de la extension cilindrica de p q (y) con x, para obtener p q (x, y). 

• Aplicar una operacion de implicacion J: J (p p (x,y), p q (x,y)) = p p _» q (x,y). 


2. Sea un sistema inteligente basado en logica difusa que maneja la calefaccion de una 
vivienda. Este sistema inteligente se encarga de calcular la temperatura que debe tener el 
agua que circula por los radiadores en funcion de la temperatura exterior y de la horn. 
Por tanto, se manejan las siguientes variables: 

• TE: Temperatura ambiente exterior, con tres valores posibles: Frio (F), Templado (T) y 
Calido (C). Su intervalo de valores es [-10,40] °C. 

• TA: Temperatura del agua que circula por los radiadores, con tres posibles valores: Agua 
Fria (AF), Agua Templada (AT) y Agua Caliente (AC). Su intervalo de valores es (0,100) 
°C. 

• TIME: horn del dia, con cuatro posibles valores: Madrugada (M), Manana (Ma), Tarde 
(T) y Noche (N). Su Intervalo de valores es [0,24). 
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Se pide: Representar graficamente las distribuciones de posibdidad de los valores 
cualitativos asociados a las variables TE y TIME. 


En este caso, se definen parcialmente las variables que debe manejar un sistema basado en el 
conocimiento para controlar el sistema de calefaccion de una vivienda, y se pide que se 
complete con la definicion de las funciones de distribucion de posibilidad asociadas a los 
valores cualitativos de dos de sus variables. Existen varias posibles soluciones, una de las 
cuales se muestra graficamente en la figura 2.12. 

^ b F ( TE ) M TE > ( TE ) Ag M (TIME)g Ma (TIME) ^ (TIME) Mnj (TIMI.) 

teQL 

-10 0 +10 +20 +30 +40 TE °C 0 3 6 9 12 15 18 21 24 TIMEh. 

Figura 2.12. Representacion grafica de las funciones de distribucion de posibilidad 
asociadas a los valores cualitativos de las variables TE y TIME. 

Es importante observar el caso de la funcion de distribucion de posibilidad asociada al valor 
cualitativo Noche (N) de la variable TIME (en trazo grueso), el cual comienza en TIME =18 
h. y termina en TIME = 3. Es decir, tiene forma de trapecio. 



3. Con respecto a las funciones de distribucion de posibilidad definidas en el ejemplo 2, ihay 
posibilidad distinta de cero de que la temperatura exterior (TE) sea a la vez fria (F) y 
tempiada (T)?. En caso afirmativo, indique para que valores de la variable TE se cumple 
esta condicion. 

SI, V TE E (0, 10). Notese que el intervalo es abierto por los dos extremos. 


4. iCual es la posibilidad de que la temperatura exterior (TE) sea templada (T) si hay 40 
°C?. i Cual es la posibilidad de una temperatura exterior (TE) calida (C) si hay 40 °C?. 

pi (40) = 0; pc (40) = 1; dado que el dominio de valores numericos para TE es cerrado. 

5. iCual es la posibilidad de que la temperatura exterior (TE) sea templada (T), teniendo en 
cuenta que hay 20 °C en el exterior y la temperatura del agua (TA) es actualmente de 40 
°C?. 

Dado que no esta definida la funcion de distribucion de posibilidad pi (TE, TA), sino solo 
para la variable TE, se calcula la extension cilindrica de pi(TE) con TA, de tal forma que: 

V TA pi (TE, TA) = pi (TE), por definicion de la extension cilindrica. Puesto que: 

Pt (TE = 20) = 1, entonces pt (TE = 20, TA = 40) = 1. De hecho, V TA pi (TE = 20, TA) = 1. 
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3.Razonamiento en logica borrosa 

Un sistema basado en el conocimiento que emplea la logica borrosa como fonnalismo 
de representacion debe definir, para un detenninado dominio, el conjunto de variables a 
emplear y su significado (semantica). Para cada una de estas variables, se define, asimismo, 
el dominio de valores numerico que puede tomar, su unidad de medida y el conjunto de 
valores cualitativos, etiquetas lingulsticas o conjuntos borrosos que se pueden asignar. Por 
ultimo, para cada valor cualitativo que puede tomar cada variable, se define su funcion de 
pertenencia o funcion de distribucion de posibilidad. 

As! mismo, se debe definir el conjunto de reglas borrosas, incluidas en la base de 
conocimiento del sistema inteligente, que relacionan a las variables entre si mediante el 
empleo de los operadores logicos borrosos: conjuncion, disyuncion, negacion o implicacion. 
Una regia borrosa sigue la siguiente sintaxis: 

antecedente —+ consecuente 

El antecedente esta formado por conjunto de clausulas de la forma variable = valor 6 
variable es valor, donde valor puede ser numerico (nltido) o cualitativo (borroso). En el 
primer caso, el valor debe pertenecer al dominio de valores numericos definido para la 
variable. En el caso de que el valor asignado a la variable sea cualitativo, se debe definir su 
correspondiente funcion de distribucion de posibilidad. Las clausulas del antecedente se unen 
mediante los operadores logicos borrosos de conjuncion (a) o disyuncion (v). Asimismo, 
cada clausula puede estar afectada por el operador de negacion (“j. El consecuente esta 
generalmente formado por una unica clausula de la forma variable = valor 6 variable es 
valor. El antecedente y el consecuente de una regia borrosa se relacionan entre si mediante el 
operador logico borroso de implicacion (—>■). Por tanto, una regia borrosa puede tener 
cualquiera de las cuatro formas siguientes, donde opi, ..., op n -i son los operadores borrosos 
de conjuncion o disyuncion, [N] representa la opcion de que la clausula siguiente pudiera 
estar afectada por la negacion y —► representa el operador borroso de implicacion: 

1. Si [N] xi es Ai opi [N] X 2 es A 2 op 2 ... op n -i [N] x n es A n , entonces y es B. 

2. [N] xi = Ai opi [N] X 2 = A 2 op 2 ... op n -i [N] x n = A n —> y = B. 

3. [N] Ai (xi) opi [N] A 2 (X 2 ) op 2 ... op n -i [N] A n (x n ) —> B (y). 

4. [N] PA 1 (xi) opi [N] PA 2 (X 2 ) OP 2 ... opn-l [N] pAn (x n ) —* pii (y). 

La primera de las formas es la mas cercana al lenguaje natural. Un ejemplo es el 
siguiente: 

Si la cantidad de CO 2 es elevada y la tendencia es positiva 0 la visibilidad no es clara, 
entonces la potencia de ventilacion es alta. 

La segunda se fonnaliza mas, empleando para ello unicamente notacion logico- 
matematica, de la fonna variable = valor (nltido 0 cualitativo): 
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(c antidad de CO 2 = elevada) a ( tendencia = positiva) v (-■ visibilidad = clara) —> 
potencia de ventilation = aha. 

La tercera forma de expresar una regia borrosa se emplea generalmente en 
controladores borrosos. En este caso, se expresa expllcitamente que los valores cualitativos se 
representan mediante funciones (de distribucion de posibilidad) que dependen de su 
correspondiente variable: 

elevada (cantidad de CO2) a positiva (tendencia) v --clara (visibilidad) —* alta (potencia 
de ventilacion). 

Por ultimo, la misma regia expresada mediante las funciones de distribucion posibilidad 
que definen los valores cualitativos (0 nitidos) que toman las variables (opcion 4): 

Peievada (cantidad de CO2) a p pos itiva (tendencia) v -qidara (visibilidad) —> p a ita (potencia de 
ventilacion). 

Para cualquier numero de clausulas que compongan el antecedente de una regia 
borrosa, siempre es posible obtener la funcion de distribucion de posibilidad, p P (x), con x 
vector, resultado de aplicar las operaciones t-nonna T (conjuncion), t-cononna S (disyuncion) 
y negacion N a las funciones de distribucion de posibilidad asociadas a los valores 
cualitativos de las clausulas del antecedente de la regia. Si fuera necesario, se deberian 
calcular previamente las correspondientes extensiones cilindricas con el fin de operar 
distribuciones de posibilidad dependientes de las mismas variables. Teniendo esto en cuenta, 
no se pierde generalidad al expresar una regia borrosa como: pa (x) —> pb (y), 0 en cualquiera 
de sus tres formas restantes. 

Asi pues, dada una base de conocimiento general fonnada por reglas borrosas de la 
forma pa (x) —> pb (y), el proceso de razonamiento 0 inferencia en logica borrosa consiste en 
calcular la funcion de distribucion de posibilidad pB’ (y) a partir de los hechos conocidos, 
expresados como pa (x). 

La logica proposicional proporciona un mecanismo de inferencia en este sentido, el 
modus ponens, que se expresa como: 

Regia: A —> B 
Hecho: A 


Resultado: B 

Dada la regia cierta A —» B, si se puede afirmar que se cumple A (el hecho conocido), 
entonces se deduce que se cumple B. 

En logica borrosa, se puede generalizar el modus ponens, al no ser necesario que se 
cumpla exactamente A, sino que es suficiente con que haya una interseccion no vacia entre el 
antecedente de la regia A y el hecho que se cumple A’, dado que ambos son conjuntos 
(borrosos). El modus ponens generalizado se podria notar como: 
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Regia: A —> B 
Hecho: A’ 


Resultado: B’ 

Teniendo en cuenta que ahora (logica borrosa) A, B, A’ y B’ son conjuntos borrosos, 
dada la regia A —» B, si no se cumple exactamente A, sino un hecho “parecido” A’ (con 
interseccion no vacla con A), no tiene porque cumplirse exactamente B, sino una afinnacion 
“parecida”, que se denomina B’. Empleando las funciones de distribucion de posibilidad que 
definen los correspondientes conjuntos borrosos A, B, A’ y B’, el modus ponens generalizado 
se puede expresar como: 

Regia borrosa: pa (x) —> p» (y) 

Hecho: pA’ (x) 


Resultado: pir (y) 

Dado que los conjuntos borrosos o etiquetas linguisticas se definen mediante sus 
correspondientes funciones de distribucion de posibilidad, es necesario, por tanto, un 
procedimiento que permita calcular la funcion de distribucion de posibilidad resultado del 
proceso de inferencia, pB’ (y), a partir de una base de conocimiento de reglas de la forma 
Pa (x) —> pb (y), conociendo que se cumple el hecho A’, representado por p^ (x). Este 
procedimiento se denomina Regia Composicional de Inferencia (RCI). Se muestra a 
continuacion como se obtiene la expresion de la regia composicional de inferencia para el 
caso de que la base de conocimiento este formada por una unica regia y se generaliza, 
posteriormente, al caso de bases de conocimiento con n reglas. 

Sean X e Y dos conjuntos nitidos (clasicos) y una relacion entre ambos, delinida por la 
funcion/(x) tal que V x G X se puede inferir y GY porque/ (x) = y: 

/: X —> Y 
x ->/(x) 

Entonces, dado un subconjunto AGX, se puede inferir, mediante la funcion/ (x), BGY 
como: B = { y GY / y =/(x), x GA}. La figura 3.1 muestra graficamente este proceso de 
inferencia: 



Figura 3.1. Inferencia del conjunto B C Y a partir de/(x) y A C X. 
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Sea R la relacion entre los conjuntos X e Y dcfinida en la llgura 3.2, donde xGX esta 
relacionado con yGY mediante la relacion R si, y solo si, <x,y> GR. Dado el subconjunto 
AGX, se puede inferir el conjunto BCY mediante la relacion R corno B = { yGY / <x,y>GR, 
xGA}. La llgura 3.2 muestra graficamente este proceso de inferencia. 



Figura 3.2. Inferencia del conjunto B C Y a partir de R y A C X. 

Se define la funcion caracterlstica de un conjunto ACX, y se nota por Pa(x), a la siguiente 
funcion de pertenencia: 

Pa: A —> {0,1} 
x —> Pa(x) 

I ^ sixeA 

P a( x H n 

0, en otro caso 

La funcion caracteristica Pa(x) devuelve 1 si x pertenece al conjunto A, y devuelve 0 en 
otro caso. Por tanto, se puede pensar que es la version nitida de la funcion de distribucion de 
posibilidad | 1 a (x). Sean A, B, la relacion R y sus respectivas funciones caracteristicas Pa(x) y 
PR(x,y) definidas en la figura 3.3.a. El procedimiento para calcular PB(y) a partir de Pa(x) y 
PR(x,y) es el siguiente, el cual se representa graficamente en las figuras 3.3.b y 3.3.c: 

1. Se proyecta la funcion caracteristica Pa(x) sobre el piano (x,y). Es decir, se hace una 
extension cilindrica de Pa(x) con y, obteniendo PA(x,y) / V y PA(x,y) = Pa(x). En la figura 
3.3.b se muestra una banda vertical de unos sobre el conjunto A. 

2. Se calcula min (PA(x,y), PR(x,y)), V (x,y). Graficamente, la figura 3.3.C muestra una 
una supcrlicic con unos mas gruesos en la interseccion entre R y la banda vertical del 
conjunto A extendi do con y. 
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3. El supremo (maximo valor de un conjunto) Vx del resultado anterior (Sup x ) devuelve 
como resultado, precisamente Pb(y) = Sup x {min (Pa(x,y), P R (x,y))}, V (x,y), lo cual se 
muestra en la figura 3.3.C. 



Figura 3.3. Proceso grafico de calculo de Pb(y) a partir de Pa(x) y Pr(x,y). a) 
Representacion de Pa(x) y PR(x,y). b) Extension de Pa(x) con y, Pa(x, y). c) Min (Pa(x, y), 
PR(x,y)), que representa la interseccion, y calculo de PB(y). 

Partiendo de la expresion obtenida en el paso 3 para el calculo de Pb(y), y 
generalizando las funciones caracterlsticas (funciones de pertenencia para conjuntos nltidos) 
a funciones de distribucion de posibilidad, se tiene: 

|i B (y) = Sup x {min (ja A (x,y), p R (x,y))}, V (x,y) 

La expresion min (j.La(x,y), p R (x,y)) proviene del calculo de la interseccion entre R y A 
extendido con y, lo cual se puede generalizar en logica borrosa mediante el empleo de una t- 
norma, T, por tanto: 

fte(y) = Sup x {t (|i A (x,y), MXy))}, v ( x ,y) 

La relacion R puede ser una relacion de implicacion en logica borrosa en donde se hacen 
corresponder puntos del conjunto A con puntos del conjunto B: 

|I B (y) = Sup x {T ((I,\(x,y), p A ^ B (x,y))}, V (x,y) 

Dado que T representa la interseccion entre conjuntos borrosos, |lA(x,y) se puede generalizar 
como |l A (x,y), donde A y A’ tienen interseccion no vacia, es decir, 3 (x,y) / |I A a a (x,y) ^ 0. 
En ese caso, se obtiene |i B ’(y) en lugar de |i B (y). 

M'fi’(y) = Sup x {T (|i\ (x,y), p A _ B (x,y))}, V (x,y) 

Por ultimo, en lugar de |l A (x,y), que seria lo ortodoxo, se suele notar como fEv(x) a la 
distribucion de posibilidad del hecho conocido, que es como originalmente suele venir 
expresado en un problema de inferencia. Por tanto, hay que tener en cuenta que la operacion 
T ( |1a(x), |i\^is(x,y)) no es posible sin realizar previamente la extension cilindrica de jl A (x) 
con y. Obviando asimismo V (x,y), dado que las operaciones con conjuntos borrosos se 
calculan siempre para todos los puntos del dominio de sus funciones de distribucion de 
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posibilidad asociadas, la expresion para la regia composicional de inferencia, que permite 
calcular la funcion de distribucion de posibilidad pB’ (y) a partir de la regia pa (x) —> pB(y) y 
del hecho, representado por pA’ (x), es la siguiente: 

My) = Sup x {T (|d/v(x), |i A ^ B (x,y))} 
a. Regia composicional de inferencia 

A partir de la formula de la regia composicional de inferencia, se observa que es 
necesario realizar las siguientes operaciones en orden: 

1. Calculo de la implicacion J para obtener la funcion de distribucion de posibilidad 
p A ^B(x,y). Dado que la regia, en la base de conocimiento, tiene la forma (expresada 
mediante las distribuciones de posibilidad asociadas a los valores cualitativos) pa (x) —> 
p B (y), no es posible operar directamente J (pa (x), pB(y)) puesto que estas dos funciones 
de distribucion de posibilidad no estan expresadas en funcion de la misma variable. Por 
ello: 


a. Se calcula la extension cilindrica de pa (x) con y para obtener pa (x,y). 

b. Se calcula la extension cilindrica de pB (y) con x para obtener pB (x,y). 

c. Se calcula J (pa (x,y), pB (x,y)), para obtener pA^B(x,y), empleando para ello la 
funcion de implicacion J elegida para el sistema inteligente construido. 

2. Calculo de la t-norma T (pa’(x), pA-m(x,y)). De nuevo, dado que la distribucion de 
posibilidad del hecho conocido, pA’(x), no esta expresado en funcion de las mismas 
variables que pa >is(x,y), es necesario: 

a. Calcular la extension cilindrica de pa’(x) con y, para obtener pA’(x, y). 

b. Calcular T (pA’(x,y), pA^B(x,y)) = pa’a(a^b) (x,y), empleando la funcion t-norma 
definida en el sistema basado en el conocimiento construido. 

3. Calculo del Sup x {pa’a(a-»b> (x,y)} = p B ’(y). 

b. Motor de inferencia borroso: generalizacion a n reglas 

Hasta ahora, todo el desarrollo de inferencia con logica borrosa, se ha realizado 
suponiendo que la base de conocimiento consta de una unica regia, expresada mediante 
funciones de distribucion de posibilidad, de la forma pa (x) —* pB(y). Si bien la fonna de 
expresar la regia se puede considerar general, dado que siempre es posible operar (o, en su 
caso, componer) todas las clausulas existentes en el antecedente de la regia, no lo es el hecho 
de que haya una unica regia en la base de conocimiento del sistema inteligente. Por ello, se 
considera ahora un sistema inteligente basado en reglas borrosas con la siguiente base de n 
reglas generales expresadas mediante funciones de distribucion de posibilidad: 
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Ri: |iAi (x) —► (iiii (y) 

R 2 : (IA2 (x) —>■ |IB2 (y) 

R n : jlAn (x) —> (J.Bn (y) 

Se supone, ademas que el hecho conocido es A’, definido como pA’(x), y se desea averiguar 

pB-(y). Para ello: 

1 . Se aplica la regia composicional de inferencia a cada regia R;: 

pBi<y) = Sup x {T (pa’(x), p A i^Bi(x,y))} 

Obteniendose parly), pB2’(y), pBn’(y). Si para alguna regia Ri ocurre que T (pA’(x), 
|XAi^Bi(x,y)) = 0, V x, y (notese la necesidad de realizar previamente la extension 
cillndrica de pa (x) con y), se obtiene como resultado que pBi’(y) = 0, V y. Dicho de otro 
modo, la regia no se dispara. Por tanto, Ri no tiene influencia en el resultado final del 
proceso de inferencia. 

2. Se calcula la funcion de distribucion final pB’(y) mediante la aplicacion de la funcion t- 
cononna S definida en el sistema inteligente a las funciones de distribucion de posibilidad 
obtenidas en el paso anterior pBr(y), pB2’(y), pBn’(y). Se puede aplicar en cualquier 
orden, dado que la t-conorma cumple las propiedades asociativa y conmutativa. No es 
necesario calcular la extension cillndrica porque todas las pBi (y) dependen de la misma 
variable y. 

PB-(y) = S (pBi’(y), S (p B 2’(y), S (p Bn (y)))...) = (J p ; (y) 

Vi=l..n 


c. Ejemplo 

Se desea construir un sistema inteligente basado en logica borrosa que permita la 
estimacion del precio de la carrera de un taxi en funcion de la distancia recorrida. Para ello, 
inicialmente, se manejan las siguientes variables: 

• D: distancia recorrida con el taxi, con posibles valores: corta (Co), nonnal (No), alrededor 
de 30 (~30), larga (L) y muy larga (ML). Su intervalo de valores es [0,100] Km. 

• P: precio de la carrera, con posibles valores: barato (B), nonnal (N) y caro (C). Su 
intervalo de valores es [0, 50] €. 

Los valores cualitativos de estas variables tienen asociadas las funciones de distribucion 
de posibilidad mostradas en la figura 3.4. Para representar el conocimiento existente acerca 
del dominio de trabajo, se dispone actualmente de una base de conocimiento formada por una 
unica regia: 

Ri: distancia = larga —> precio = caro 

Rl: piarga (D) * Llca:o (P) 
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Figura 3.4. Funciones de distribucion de posibilidad asociadas a los valores cualitativos 

de las variables D (distancia) y P (precio). 

Se pide: suponiendo que se desea tomar un taxi para recorrer una distancia de alrededor de 30 
Km. (D = -30), definido por (j~3o (distancia), aplicar la regia composicional de inferencia 
(RCI) con la t-norma y t-conorma de Zadeh para estimar los valores numerico y cualitativo 
del precio que costara la carrera. 

SOLUCION: 

Se muestran a continuacion los operadores a utilizar, teniendo en cuenta que el 
enunciado indica que se deben utilizar la t-norma y t-conorma de Zadeh y, por tanto, una de 
las dos funciones de implicacion J que se deducen de ellas. Asi, V x, con x, y dos funciones 
de distribucion de posibilidad generales, p P (x) y p q (x) respectivamente: 

• T (x, y) = min (x, y). 

• S (x, y) = max (x, y). 

• N (x) = 1- x. 

• J (x, y) = max (1 - x, y), que es una de los funciones de implicacion que se deducen. 

Dado que el ejercicio maneja unicamente dos variables, P y D, se podria resolver de 
forma grafica, lo cual daria como resultado un resultado exacto. Sin embargo, el 
procedimiento grafico solo es valido en problemas donde, como este, se manejan dos 
variables, no siendo, por tanto, general. Por el contrario, se opta por un desarrollo tabular que, 
si bien es inexacto por tener que discretizar las funciones de distribucion de posibilidad 
(continuas) es valido para cualquier numero de variables. 

La eleccion de los intervalos de discretizacion (regulares) de las funciones de 
distribucion de posibilidad es un asunto importante. Si se elige un radio para los intervalos de 
discretizacion demasiado grande, se obtiene una solucion demasiado inexacta. Por el 
contrario, la eleccion de un radio pequeno, proporciona un resultado mas preciso, pero con un 
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mayor coste computacional. En el ejemplo propuesto, se ha elegido un radio para los 
intervalos de discretizacion de 5 Km. para el caso de la variable D, y de 10 € para P. 

La tabla 3.1 muestra el resultado de aplicar la regia composicional de inferencia al caso 
propuesto en el problema. Las dos primeras columnas muestran, respectivamente, los valores 
numericos en donde se discretizan las funciones de distribucion de posibilidad que dependen 
de las variables PyD, denominados puntos de discretizacion. Las columnas 3 y 4, pi ar ga (D,P) 
y pcaro (D,P), muestran los valores de posibilidad en cada uno de los puntos de discretizacion. 
Para poder calcular piarga-^caro (D,P), es necesario hacer las extensiones cilindricas de piarga (D) 
con P y pcaro (P) con D para obtener, respectivamente, piarga (D,P) y p ca ro (D,P). En el caso de 
piarga (D,P), la extension cilindrica se representa en la tabla mediante la repeticion de los 
valores de posibilidad en todas las celdas de la columna piarga (D,P) para cada uno de los 
valores de precio (P). Efectivamente, si la extension cilindrica significa, por definicion, que V 
P piarga (D,P) = piarga (D), es necesario copiar los valores de posibilidad de piarga (D) de cada 
punto de discretizacion de D en todos los puntos de discretizacion de P. En el caso de p ca ro 
(D,P), que solo depende de precio, bastaria con anotar en la cuarta columna, p ca ro (D,P), para 
cada punto de discretizacion de D, el unico valor de posibilidad correspondiente al punto de 
discretizacion de P. 

A partir de piarga (D,P) y p ca ro (D,P), es posible calcular piarga^caro (D,P) para cada valor 
de posibilidad de las columnas piarga (D,P) y p ca ro (D,P). El resultado se muestra en la quinta 
columna, empleando para ello J (x, y) = max (1 - x, y), V D, P, tal y como se deduce del 
enunciado. En este caso, x = piarga (D,P), e y = p ca ro (D,P). Por ejemplo, para D = 25 y P = 0, 
se tiene que piarga (25, 0) = 0,5 y p ca ro (25, 0) = 0. Entonces, J (piarga (25, 0), p ca ro (25, 0)) = 
max (1 - 0,5, 0) = 0,5, que es el resultado mostrado en la sexta fda de la primera celda de la 
quinta columna. 

La sexta columna representa los valores de posibilidad de [J~30 (D), extendido 
cilmdricamente con P, de forma similar al caso de piarga (D,P) en la tercera columna. Dado 
que piarga—>car 0 (D,P) y la extension cilindrica de (i ~30 (D) con P, mostrada en la columna sexta, 
dependen de las mismas variables D y P, es posible calcular la t-norma de ambas funciones 
mediante la expresion indicada en el enunciado: T (x, y) = min (x, y), V D, P. El resultado se 
muestra en la columna septima, llamada T. De nuevo, en el caso de la sexta fila de la primera 
celda de la columna T, se tiene: T (piarga —►caro (25,0) , [_l~30 (25,0) = min (0,5, 1) = 0,5. 

Finalmente, la ultima columna muestra los valores de posibilidad de la funcion de 
distribucion resultado pb- (P) para cada uno de los puntos de discretizacion elegidos para la 
variable P. Para ello, se calcula, para cada valor de P, el maximo de las posibilidades de la 
funcion T (columna anterior), para todo D (Supo). En el caso de la primera celda de la ultima 
columna, el valor 0,5 es el maximo valor de posibilidad de T para P = 0 y cualquier valor de 
D, es decir, SupD (T (D, 0)) = 0,5. 
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Tabla 3.1. Resultado de aplicar la regia composicional de inferencia al ejemplo 


propuesto. 


P: Precio 

D: Distancia piarga (D,P) 

jlcaro (D,P) 

j-llarga * caro (D,P) 

p-30 (D,P) 

T Supd: gB- (P) 


0 

0 

0 

i 

0 

0 : 


5 

0 

0 

i 

0 

0 


10 

0 

0 

i 

0 

0 


15 

0 

0 

i 

0 

0 


20 

0 

0 

i 

0 

0 


25 

0,5 

0 

0,5 

1 

0,5 

0 

30 

1 

0 

0 

1 

0 0,5 


35 

0,5 

0 

0,5 

1 

0,5 


40 

0 

0 

1 

0 

0 


45 

0 

0 

1 

0 

0 i 


50 

0 

0 

1 

0 

o : 


100 

0 

0 

1 

0 

0 


0 

0 

0 

1 

0 

0 


5 

0 

0 

1 

0 

0 


10 

0 

0 

1 

0 

0 


15 

0 

0 

1 

0 

0 


20 

0 

0 

1 

0 

0 


25 

0,5 

0 

0,5 

1 

0,5 

10 

30 

1 

0 

0 

1 

0 0,5 


35 

0,5 

0 

0,5 

1 

0,5 


40 

0 

0 

1 

0 

0 


45 

0 

0 

1 

0 

0 


50 

0 

0 

1 

0 

0 


100 

0 

0 

1 

0 

0 


0 

0 

0 

1 

0 

0 


5 

0 

0 

1 

0 

0 


10 

0 

0 

1 

0 

0 


15 

0 

0 

1 

0 

0 i 


20 

0 

0 

1 

0 

0 


25 

0,5 

0 

0,5 

1 

0,5 

20 

30 

1 

0 

0 

1 

0 0,5 


35 

0,5 

0 

0,5 

1 

0,5 


40 

0 

0 

1 

0 

0 


45 

0 

0 

1 

0 

0 


50 

0 

0 

1 

0 

0 


100 

0 

0 

1 

0 

0 


0 

0 

0,5 

1 

0 

0 


5 

0 

0,5 

1 

0 

0 


10 

0 

0,5 

1 

0 

0 


15 

0 

0,5 

1 

0 

0 


20 

0 

0,5 

1 

0 

0 


25 

0,5 

0,5 

0,5 

1 

0,5 

30 

30 

1 

0,5 

0,5 

1 

0,5 0,5 


35 

0,5 

0,5 

0,5 

1 

0,5 


40 

0 

0,5 

1 

0 

0 


45 

0 

0,5 

1 

0 

0 


50 

0 

0,5 

1 

0 

0 


100 

0 

0,5 

1 

0 

0 
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P: Precio 

D: Distancia 

piarga (D,P) 

Pcaro (D,P) glarga—> caro (D,P) 

p-30 (D,P) 

T 

SupDl PB’ (P) 


0 

0 

1 

1 

0 

0 


5 

0 

1 

1 

0 

0 


10 

0 

1 

1 

0 

0 


15 

0 

1 

1 

0 

0 


20 

0 

1 

1 

0 

0 


25 

0,5 

1 

1 

1 

1 

40 

30 

1 

1 

1 

1 

1 

1 


35 

0,5 

1 

1 

1 

1 


40 

0 

1 

1 

0 

0 


45 

0 

1 

1 

0 

0 


50 

0 

1 

1 

0 

0 


100 

0 

1 

1 

0 

0 


0 

0 

1 

1 

0 

0 



5 

0 

1 

1 

0 

0 



10 

0 

1 

1 

0 

0 



15 

0 

1 

1 

0 

0 



20 

0 

1 

1 

0 

0 



30 

0,5 

1 

1 

1 

1 


50 

35 

1 

1 

1 

1 

1 

1 


40 

0,5 

1 

1 

1 

1 



45 

0 

1 

1 

0 

0 



50 

0 

1 

1 

0 

0 




0 

1 

1 

0 

0 



100 









0 

1 

1 

0 

0 



La figura 3.5a. muestra graficamente los valores de posibilidad de pir (P) en los puntos de 
discretizacion de P, obtenidos como resultado de aplicar la regia composicional de inferencia, 


tal y como se observa en la ultima columna de la tabla 3.1. La figura 3.5.b muestra 
graficamente pir (P) despues de unir los valores de posibilidad mediante rectas (interpolacion 
mediante polinomios de primer grado). Notese que cuanto mayor sea el numero de puntos de 
discretizacion, mas precisa es la funcion pir (P) que se obtiene como resultado, aunque como 
consecuencia de un mayor numero de calculos. 


A 

1 a) * 


0,5 


* 


* 


* 


0 10 20 30 40 


* 


Precio 

-> 

50 € 



Figura 3.5. Representacion grafica de pB’ (P), resultado de aplicar la RCI. 


d. Proceso de interpretation del resultado 

El resultado obtenido en el ejercicio anterior, pir (P), debe ser interpretado. Se debe 
saber cual es la semantica del resultado, cuanto costara el taxi si se recorre una distancia de 
alrededor de 30 Km. (D = ~30). De hecho, en el ejercicio se pide que se estimen los valores 
numerico y cualitativo del precio que costara la carrera. Por tanto, es necesario convertir la 
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funcion de distribucion de posibilidad obtenida como resultado del proceso de inferencia en 
un valor numerico o un valor cualitativo de entre los conocidos: barato (B), normal (N) o caro 
(C). Este proceso de conversion, ya sea en un valor numerico o cualitativo, se denomina en 
ingles defuzzify , que se podrla espanolizar usando el barbarismo desborrocificacion, que no 
existe. Otra opcion es realizar una traduccion libre del termino en base a su semantica en el 
dominio de la logica borrosa: interpretar o convertir. 

Existen diversos procedimientos para convertir una funcion de distribucion de 
posibilidad en un valor numerico o cualitativo. Para el primer caso, se describe el metodo del 
centro de gravedad gir. Para el calculo de valores cualitativos, se exponen los metodos de la 
distancia y ^cuenta, leldo como “sigma cuenta”. 

El metodo del centro de gravedad convierte una funcion de distribucion de 
posibilidad en un valor numerico a partir, precisamente, del calculo de su centro de gravedad. 
Dada una funcion de distribucion de posibilidad generica p p (x) y un conjunto de puntos de 
discretizacion de la funcion {xi, X 2 , ..., x n }, elegidos, recuerde, de forma que sean 
equidistantes, el centro de gravedad de p P (x) se define como: 

n 

XvPpOo 

g P = - 

Xm x .) 

i-l 

En el caso del ejemplo de la figura 3.5, el centro de gravedad de pir (P), que permite 
calcular el valor numerico del precio de una carrera de alrededor de 30 Km. de distancia, se 
calcula de la siguiente forma, teniendo en cuenta que se toma un punto de discretizacion cada 
10 €. 

0-0,5 + 10*0,5 + 20-0,5 + 30-0,5 + 40-1 + 50.1 120 

g R1 =--=-= 30€ 

B 0,5 + 0,5 + 0,5 + 0,5 + 1 + 1 4 

Es decir, si se recorre una distancia de alrededor de 30 Km., el taxi costara 30 €. Si esta 
prediccion coincide con la realidad, el sistema basado en el conocimiento esta bien 
construido para este caso. Por el contrario, si el resultado es erroneo, es necesario corregir el 
sistema inteligente modificando la base de reglas, las funciones T, S, N y J, las funciones de 
distribucion de posibilidad o el metodo de conversion a valores numericos. La figura 3.6 
muestra graficamente el centro de gravedad de pir (P). 
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Es importante tener en cuenta que, tal y como se define en el ejercicio propuesto en la 
seccion 3.c , los valores numericos de la variable precio E [0,50] €. Es decir, tanto el 0 como 
el 50 pertenecen al intervalo, son valores que la variable puede tomar. Por esa razon, estos 
dos valores deben ser incluidos en el calculo del centro de gravedad de pB’ (P). No solo eso, 
sino que, ademas, su inclusion influye en el resultado final: si bien es cierto que el valor 0 no 
influye en el resultado del numerador para el calculo del centro de gravedad, si Xi = 0, Xi • p p 
(xi) = 0, no ocurre lo mismo en el caso del denominador, en el ejemplo, p p (x;) = 0,5 para Xi = 
0 . 


El metodo de la distancia devuelve el valor cualitativo Ai, de entre los conocidos o 
definidos en el dominio de valores cualitativos de la variable en cuestion (en el ejemplo: 
barato, normal o caro), cuya funcion de distribucion de posibilidad se parece mas a la funcion 
de distribucion de posibilidad obtenida como consecuencia de un proceso de inferencia en 
logica borrosa, notada como pir (x). Para ello, se define una distancia entre dos funciones de 
distribucion de posibilidad que mide la similitud o grado de parecido entre ellas. Este metodo 
consiste en calcular la distancia de pB’ (x) a cada una de las funciones de distribucion de 
posibilidad definidas para cada valor cualitativo de la variable x, devolviendo como resultado 
aquel Ai cuya funcion de distribucion de posibilidad minimiza d (pir (x), pAi (x)) V i. 

^Cuando podemos decir que dos distribuciones de posibilidad, pB' (x) y pAi (x), son 
parecidas?. Una posible respuesta es: cuando sus correspondientes centros de gravedad sean 
cercanos. Efectivamente, las figuras 3.7.a y 3.7.b muestran, respectivamente, el caso de dos 
distribuciones de posibilidad parecidas (con centros de gravedad similares) y dos 
distribuciones de posibilidad que no lo son. Sin embargo, esta medida no es suficiente como 
se puede observar en la figura 3.7.C donde las dos distribuciones de posibilidad tienen el 
mismo centro de gravedad, pero no se pueden considerar que sean parecidas. 



Figura 3.7.a) dos funciones de distribucion de posibilidad parecidas con centros de 
gravedad similares. b) dos funciones de distribucion de posibilidad que no son 
parecidas, al igual que sus centros de gravedad. c) dos funciones de distribucion de 
posibilidad que no siendo parecidas, tienen el mismo centro de gravedad. 

A partir de la figura 3.7, sobre todo en el caso c), se observa que el calculo de los 
centros de gravedad no es una medida suficiente para medir el grado de parecido o distancia 
entre dos funciones de distribucion de posibilidad. Se podrla pensar que las dos 
distribuciones de la figura 3.7.C no son parecidas porque la medida de sus areas son muy 
dispares. Pero la medida del area de una distribucion de posibilidad tampoco es suficiente, 
como se puede observar en el caso de la figura 3.7.b, donde las dos distribuciones tienen el 
mismo area, pero no son parecidas. 
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Se puede entonces pensar en una medida de distancia entre dos funciones de 
distribucion de posibilidad d (pir (x), pAi (x)) que combine tanto su centra de gravedad como 
su area. En este caso, dos funciones de distribucion de posibilidad son cercanas entre si (o 
parecidas ), estan a corta distancia, si tanto sus centros de gravedad como sus areas son 
similares. Efectivamente, en este caso, las distribuciones de la figura 3.7.a son cercanas y, por 
tanto, parecidas, mientras que las funciones de las figuras 3.7.b y 3.7.C no lo son. Una 
formula para la distancia d (pir (x), pAi (x)) que tenga en cuenta estos dos aspectos puede ser 
la siguiente: 

d(M*),Mx)) ~ V a ‘ (g B' _ ^ A i P (^B 1 _ ) 

Donde Ab- y Aa. representan, respectivamente, el area de las funciones de distribucion 
de posibilidad pB’ (x) y pAi (x), correspondientes a B’ y Ai. gB’ y gAi son, respectivamente, sus 
centros de gravedad. Por ultimo, a y p son dos factores de ponderacion de la influencia que 

tienen los centros de gravedad y areas de las distribuciones de posibilidad en el resultado 
final del calculo de la distancia entre ambas: a, /? E[0,1], a + /? = 1. 

En el caso del ejemplo de la seccion 3.c, se deben calcular las distancias de la funcion 
de distribucion de posibilidad obtenida como resultado de la RCI a cada una de las tres 
distribuciones de posibilidad dclinidas para los valores cualitativos incluidos en el dominio 
de la variable precio (P): barato (B), normal (N) y caro (C). Se obtiene como resultado el 
valor cualitativo cuya funcion de distribucion de posibilidad minimiza la distancia a pB’ (P). 
Para ello, se calculan, en primer lugar, los centros de gravedad de cada una de las cuatro 
distribuciones de posibilidad (tomando un punto de discretizacion cada 10 €). 

1. gB’ = 30, calculado anteriormente. 

2. g B = (0 1 + 10-1 + 20 0) / (1 + 1 + 0) = 5. 

3. g N = (0-0+ 10-1 +20-1 +30-1 +40-0)/(0+ 1 + 1 + 1 + 0) = 20. 

4. gc = (20-0 + 30-0,5 + 40-1 + 50-1) / (0 + 0,5 + 1 + 1) = 42. 

A continuacion, se calculan las areas de cada una de las cuatro distribuciones de 
posibilidad. Para ello, se particiona la superficie de cada funcion distribucion de posibilidad 
en figuras geometricas cuyo area sea sencilla de calcular. Por ejemplo, triangulos o 
cuadrados. El area de la funcion de distribucion de posibilidad es el sumario de las areas de 
las figuras geometricas que conforman la particion. Se muestra, en la figura 3.8, a modo de 
ejemplo, una particion de la superficie de pb- (P), que consta de dos rectangulos y un 
triangulo cuyas areas son Abu = 25, Abu = 5 y Abu = 2,5 respectivamente. Por tanto, 
A B ’=32,5. 

1. Ab’ = 25 + 5 + 2,5 = 32,5 (Abu + Abu + Abu). 

2. Ab =10 + 5 15 (Ab rectangulo + Ab triangulo)- 

3. An = 5 + 30 + 5 = 40. (An triangulo 1 + An rectangulo + An triangulo2). 

4. Ac = 10 + 10 = 20 (Ac triangulo + Ac rectangulo). 
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Figura 3.8. Calculo del area de una funcion de distribucion de posibilidad. 

Por ultimo, se calculan las distancias de pB’ (P) a cada una de las tres distribuciones de 
posibilidad conocidas, tomando, por ejemplo, a = 0,6 y /? = 0,4 para dar mayor peso 

(importancia) al centro de gravedad frente al area en el calculo de la medida de parecido 
(distancia) entre dos funciones de distribucion de posibilidad: 

d (p»- (P), p B (P)) = 0,6-(30-5) 2 + 0.4-(32,5-15) 2 = 375 + 122,5 = 497,5 

d (p B ’ (P), p N (P)) = 0,6 (30-20) 2 + 0.4 (32,5-40) 2 = 60 + 22,5 = 82,5 

d (p B * (P), pc (P)) = 0,6-(30-42) 2 + 0.4-(32,5-20) 2 = 86,4 + 62,5 = 148,9 

Se ha omitido el calculo de las raices cuadradas, dado que no es necesario calcular 
exactamente el valor de cada una de las distancias, sino unicamente valores aproximados a 
efectos de comparar unas distancias con otras para elegir la menor. En este caso, la distancia 
menor es d (pB’(P), Pn(P)), es decir, pn(P) esta mas cerca (se parece mas) a pb (P). Por tanto, 
la respuesta que debe dar el sistema basado en el conocimiento, en terminos cualitativos, a 
la pregunta planteada en el ejemplo del apartado 3.c es que el precio del taxi sera normal, 
que es el valor incluido en el dominio de valores cualitativos de la variable precio (P) cuya 
funcion de distribucion de posibilidad es pn(P). 

El uso del metodo de la distancia para calcular a que funcion de distribucion de 
posibilidad (valor cualitativo) de las conocidas se parece mas la distribucion de posibilidad 
obtenida como resultado de un proceso de inferencia en logica borrosa, notado como pB’(x), 
tiene de un coste elevado en terminos de carga computacional o numero de calculos a 
realizar. Aun obviando el calculo de la raiz cuadrada, es necesario determinar cuales son las 
areas de todas las funciones de distribucion de posibilidad implicadas, sus centros de 
gravedad y calcular la suma de las diferencias al cuadrado. 

El metodo de la Xcuenta (leido como sigma-cuenta ) tiene el mismo objetivo que el 
metodo de la distancia: devolver el valor cualitativo Ai, de entre los definidos en el dominio 
de la variable en cuestion, cuya funcion de distribucion de posibilidad se parezca mas a la 
distribucion de posibilidad obtenida como consecuencia de un proceso de inferencia en logica 
borrosa, notado como pB’ (x). Sin embargo, requiere de un menor coste computacional 
porque esta basado en el calculo de estimaciones de las areas de las funciones de distribucion 
y cercania de los centros de gravedad. Por tanto, la £cuenta emplea los mismo criterios que el 
metodo de la distancia para establecer el grado de parecido o similitud entre dos funciones de 
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distribucion de posibilidad, area y centro de gravedad, pero con menor coste computacional 
porque calcula estimaciones y no valores exactos. 

Se define la £ CU enta de una funcion de distribucion de posibilidad pa (x), y se nota como 

^jcuenta(A)! 

n 

y (A)=Y[i A (x,) 

cuenta r A 1 

i=l 

Donde {xi, X 2 , x n } es el conjunto de puntos de discretizacion de la funcion, elegidos 
de forma que sean equidistantes entre si. La Xcuenta(A) se p ue de considerar una estimacion del 
area de la funcion de distribucion de posibilidad pa (x), dado que suma las alturas (valores de 
posibilidad) de la funcion de distribucion de posibilidad pa (x) en cada uno de los puntos de 
discretizacion x,. El termino estimacion no viene dado porque £ C uenta(A) ~ a a, sino porque, 
dadas dos funciones de distribucion de posibilidad pai(x), pA 2 (x) / Aai < Aa2 —+ X cuenta(Al) ~ 

^cuenta(A2). 

Tomando como ejemplo las funciones de distribucion de posibilidad pb (P), pn (P), 
pc (P) y pb’ (P) que representan los diferentes valores cualitativos de precios de 
desplazamientos en taxi del ejemplo planteado en la seccion 3.c, se tiene que los valores 
exactos para las areas, ordenadas de menor a mayor son (calculadas en la resolucion del 
ejercicio con el metodo de la distancia): Ab = 15; Ac = 20; Ab’ = 32,5; An = 40. Se calculan 
a continuacion, a modo de ejemplo, la Xcuenta de cada una de estas funciones para observar 
como, efectivamente, se mantiene, aproximadamente, el mismo orden (se ha tornado un 
punto de discretizacion cada 10 €): 

Xcuenta(B) =1 + 1+ 0 = 2; X cuenta (C) = 0 + 0,5 + 1 +1 = 2,5. 

Xcuenta(B’) = 0,5 + 0,5 + 0,5 + 0,5 + 1 + 1 = 4; £cuenta(N) = 0+l + l + l+0 = 3. 

La distancia entre los centros de gravedad de dos funciones de distribucion de 
posibilidad pai(x) y pA 2 (x) se podria estimar, inversamente, como la Xcuenta de la interseccion 
de sus valores cualitativos asociados (o conjuntos borrosos). Es decir, Xc L| ema(Ai a A 2 ). En el 
ejemplo de la figura 3.7, las funciones de distribucion de posibilidad de los casos a) y c) 
tienen sus centros de gravedad cercanos (su distancia es practicamente nula 0 cero). En estos 
dos casos, la £cuenta c i e ] a interseccion de sus valores cualitativos asociados, que se puede 
interpretar como el area de la funcion de distribucion de posibilidad resultado de la 
interseccion, es elevada, por lo que se asume que la distancia entre sus centros de gravedad es 
pequena. Por el contrario, la distancia entre los centros de gravedad de las dos funciones de 
distribucion de posibilidad de la figura 3.7.b es elevada. En este caso, la Amenta de la 
interseccion de sus valores cualitativos asociados es cero porque dicha interseccion es vacia. 

Segun el razonamiento anterior, el valor cualitativo A; cuya funcion de distribucion de 
posibilidad se parece mas a la funcion de distribucion de posibilidad obtenida como 
consecuencia de un proceso de inferencia en logica borrosa (RCI) pb- (x), deberia ser aquel 
que maximice la expresion Xcuenta(B’AAi), V Ai. Sin embargo, en los ejemplos de la figura 
3.7, donde se supone que las funciones de distribucion de posibilidad del caso a) se parecen 
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mas que las del caso c), podrla ocurrir que Xcuenta(B’ aAroJo izda.) < X^'cmalB’ aAr OJ o Doha.). Es 
decir, que la £ C uenta de la interseccion de los valores cualitativos asociados sea mayor en el 
caso de las funciones de distribucion de posibilidad del ejemplo de la figura 3.7.C que en el 
caso de la figura 3.7.a. De hecho, esto es lo que parece ocurrir a simple vista. 


Este error se produce porque AARojojzda.«AARojo Dcha.. Cuanto mayor es el area de una 
funcion de distribucion de posibilidad, mayor oportunidad tiene su correspondiente valor 
cualitativo o conjunto borroso de intersecar con cualquier otra. El area AR 0 j 0 _Dcha. en la figura 
3.7.C es tan grande que practicamente abarca todo el dominio de x. En esas condiciones, es 
logico pensar que el area de la funcion de distribucion de posibilidad resultado de la 
interseccion de RojoDcha. con B’ sea elevada. Por tanto, para calcular el grado de parecido 
entre dos funciones de distribucion de posibilidad pb- (x) y pAi (x), no solo basta con calcular 
Xcuenta(B’AAi), sino que tambien es necesario tener en cuenta, de fonna inversamente 
proporcional, el area de la funcion de distribucion de posibilidad de Ai, o bien, una 
estimacion de esta, es decir, Xcuenta(Ai). Como conclusion, dos funciones de distribucion de 
posibilidad, pb- (x) y pAi (x), son parecidas cuanto mas elevado sea el area de la interseccion 
de sus correspondientes valores cualitativos (conjuntos borrosos) A; y B’, y mas pequena el 
area de la funcion de distribucion de posibilidad con la que se compara, pAi (x). Es decir, 
Pb- (x) se parece mas a aquella pAi (x) V Ai que maximice Xcuenta (Ai / B’): 

, Y (B' A Aj) 

Y A i /B , =%aE r ^ 

cuenta ' 1 ' X ( A ) 

cuenta ' 1 ' 

Siendo pAi (x) cada una de las funciones de distribucion de posibilidad definidas para 
los valores cualitativos del dominio de las variables manejadas en el sistema inteligente 
basado en logica borrosa con el que se trabaja. 

En el caso del ejercicio del calculo del precio de una carrera de taxi, seccion 3.c, es 
necesario calcular, en primer lugar la ^cuenta de las funciones de distribucion de posibilidad de 
los valores cualitativos conocidos (barato, normal y caro): £ cuenta(B), ycuenta(N) y ^ cuenta(C)j 
lo cual ya se ha llevado a cabo tomando un punto de discretizacion cada 10 €: 


Zcuenta(B) = 2; Xcuenta(N) = 3 ; Icuenta(C) = 2,5 

Posterionnente, se calcula la ^cuenta de la interseccion de estos valores cualitativos 
conocidos con B’, obtenido como consecuencia del proceso de inferencia llevado a cabo con 
la regia composicional de inferencia: y C uenta(B’ aB), Xcuenta(B’ aN) y Xcuenta(B’ aC). La figura 
3.9 muestra graficamente las funciones de distribucion de posibilidad resultado de las 
intersecciones, utilizando la t-norma de Zadeh, T (x, y) = min (x, y), V P. 


De acuerdo con la figura 3.9 y tornado puntos de discretizacion, de nuevo, cada 10 €: 
Xcuenta(B’ aB) = 0,5 + 0,5 = 1 
Icuenta(B’AN) = 0,5 + 0,5 + 0,5 = 1,5 
Icuenta(B’AC) = 0,5 + 1 + 1=2,5 


Por ultimo, se calcula ^cuenta (Ai / B’), Ai= {B, N, C}: 
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Icuenta (B/ B’) =1/2; ^cuenta (N/ B’) = 1,5/3; X— (C / B’) = 2,5/2,5 



Figura 3.9. Calculo de B’a B, B’a N y B’a C mediante la t-norma de Zadeh. 


La funcion de distribucion de posibilidad que maximiza Xcuenta (A; / B’) es pc (P), por 
tanto, A; = C. Es decir, aquel valor cualitativo Ai que maximiza el area de la funcion de 
distribucion de posibilidad resultado de la interseccion entre Ai y B’, independientemente de 
su propia area Am. Por tanto, se considera que la funcion de distribucion de posibilidad de C 
(caro) es la que mas se parece a la distribucion de B’, resultado de aplicar el proceso de 
inferencia (regia composicional de inferencia). Ante la pregunta planteada en el ejemplo del 
apartado 3.c: <;,cual sera el precio (cualitativo) de la camera de un taxi si la distancia a recorrer 
es de alrededor de 30 Km.?; la respuesta es, de acuerdo con el metodo de la Xcuenta, que el 
taxi sera caro. 

Se han aplicado dos metodos diferentes para el calculo del valor cualitativo del precio 
de la camera de un taxi, suponiendo que se desean recorrer alrededor de 30 Km, y se han 
obtenido dos resultados diferentes. ^Cual es entonces el resultado correcto?. La respuesta es, 
en terminos formales, que los dos son correctos. Seria necesario validar estos resultados en el 
entomo real de trabajo para averiguar cual de las dos soluciones es mas satisfactoria. 
Tambien es posible modificar los valores de ponderacion a y p de la formula del calculo de 
distancias entre dos funciones de distribucion de posibilidad. En ese caso, tambien se 
obtendrian, posiblemente, valores cualitativos diferentes, pero todos ellos formalmente 
correctos. Ademas de estos parametros, existen otras decisiones que se deben tomar a la hora 
de disenar un sistema inteligente basado en logica borrosa, como por ejemplo: la distancia 
entre los puntos de discretizacion, los operadores t-norma, t-conorma, negacion e implicacion 
(siempre que se satisfaga la condicion de dualidad), o la forma que adoptan las funciones de 
distribucion de posibilidad que definen los valores cualitativos o conjuntos borroso. Todas 
estas decisiones influyen en el proceso de inferencia, pudiendo proporcionar resultados 
diferentes para cada una de ellas, aunque formalmente correctos. El ajuste de estos 
parametros ( tunning , en ingles) se puede llevar a cabo mediante un proceso de validacion del 
sistema inteligente en el entomo real de trabajo. 
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4.Controladores borrosos 

Los controladores borrosos, o difusos, son sistemas inteligentes basados en la 
representation del conocimiento con imprecision mediante la logica borrosa que permiten 
mantener a un sistema en equilibro a pesar de las perturbaciones que pueda sufrir. Un 
controlador borroso emplea una base de conocimiento formada por reglas borrosas para 
describir el sistema a controlar, a diferencia de los controladores convencionales, que utilizan 
ecuaciones diferenciales. Los controladores borrosos reciben como entrada un conjunto de 
variables de estado que describen la situacion actual del sistema a controlar y devuelve como 
salida las acciones a llevar a cabo con el fin de mantener al sistema en su situacion de 
equilibro. La figura 4.1 muestra de forma esquematica la interaccion del controlador borroso 
con el sistema a controlar. 



Figura 4.1. Interaccion del controlador borroso con el sistema a controlar. 

Los controladores borrosos se pueden emplear en multiples aplicaciones de control o de 
mantenimiento de sistemas en equilibro en donde el tiempo de respuesta no es critico, siendo 
la facilidad de implementacion su principal ventaja. Dos ejemplos de aplicacion de los 
controladores borrosos son los siguientes: 

1. Sistemas de climatizacion: en este caso, se desea mantener un recinto, como una 
vivienda u oficina, a una temperatura constante a pesar de la temperatura exterior 
(perturbaciones extemas al sistema), dado que los sistemas de aislamiento no son 
perfectos. Si la temperatura exterior es muy elevada, la temperatura del recinto tiende a 
aumentar, por lo que es necesario que el sistema de climatizacion produzca frio para 
mantener la temperatura interior constante. Por el contrario, si la temperatura exterior es 
muy baja, el sistema de climatizacion debe producir calor. Por tanto, el controlador 
borroso se encargaria de gobemar el sistema de climatizacion. Las variables de estado que 
llegan al controlador borroso podrian ser en este caso: la temperatura actual del recinto, 
medida mediante sensores, el desvio o error con respecto a la temperatura de referenda o 
de control y la variacion o tendencia de la temperatura del recinto (aumenta, disminuye o 
es constante). La respuesta del controlador borroso, su accion, puede venir dada en forma 
de cambios en el estado del elemento a controlar, como el caudal de gas o niveles de 
apertura de valvulas. 

2. Sistemas de ventilacion: suponga que se instala un sistema de ventilacion forzada en 
un recinto cerrado, como puede ser un tunel o aparcamiento subterraneo, para mantener el 
aire limpio a pesar de las perturbaciones que producen el transito de vehiculos a motor. 
En este caso, el controlador borroso debe gobemar la potencia de ventilacion o extraccion 
de aire (accion) en funcion de (variables de estado) ciertas medidas de calidad del aire en 
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un momenta dado, desvio con respecto a lo que se considera un aire limpio y la tendencia 
de la calidad del aire. 

Los valores de las variables de estado y acciones a tomar pueden ser etiquetas 
lingulsticas o conjuntos borrosos definidos mediante sus correspondientes funciones de 
distribucion de posibilidad o funciones de pertenencia. Por tanto, en un controlador borroso 
es necesario definir, de la misma fonna que en un sistema inteligente basado en logica 
borrosa general: el conjunto de variables que se manejan, sus dominios de valores numericos, 
dominios de valores cualitativos y sus correspondientes funciones de distribucion de 
posibilidad. As! mismo, la base de conocimiento incluye un conjunto de reglas borrosas, la 
cuales definen, en sus antecedentes, posibles situaciones que se pueden producir en el sistema 
a controlar; y en los consecuentes, para cada situacion descrita, la accion a llevar a cabo para 
devolver al sistema a controlar a su situacion de equilibrio. 

Sea s una variable que mide un detenninado parametro del estado de un sistema a 
controlar (por ejemplo, la temperatura ambiente en el ejemplo del sistema de climatizacion), 
e es el error o desvio con respecto a la situacion de equilibrio, Ae la variacion o tendencia del 
error y v describe la accion a tomar (por ejemplo, calentar agua a traves de una caldera para 
mantener la temperatura ambiente). Suponiendo que se han dclinido las funciones de 
distribucion de posibilidad de los valores cualitativos que pueden tomar cada una de estas 
variables, una base de n reglas de un controlador borroso podria ser la siguiente: 

Ri: s = Ai a e = Bi a Ae = Ci —> v = Di 

Rn s = Ai a e = Bi a Ae = Ci —» v = Di 

Rn. S An A e Bn A Ae Cn ^ V Dn 

Donde A y —» representan, respectivamente, las operaciones de conjuncion 
(interseccion) e implicacion entre etiquetas linguisticas o conjuntos borrosos. Notese como 
las clausulas de los antecedentes de las reglas de una base de conocimiento de un controlador 
borroso estan siempre unidas mediante conjunciones, lo cual pennite simplificar el proceso 
de inferencia. Cada etiqueta linguistica o conjunto borroso en la base de conocimiento, Ai, Bi, 
Ci y Di, Vi, esta definida mediante su correspondiente funcion de distribucion de posibilidad, 
respectivamente, pAi(s), pBi(e), pci(Ae) y pDi(v), Vi. En controladores borrosos es comun 
emplear una notacion diferente, aunque equivalente, para expresar las funciones de 
distribucion de posibilidad. En este caso, se utiliza el nombre de la etiqueta linguistica o 
conjunto borroso para notar la correspondiente funcion de funcion de distribucion de 
posibilidad. Asi, la funcion de distribucion de posibilidad pa(x), que define a la etiqueta o 
conjunto borroso A, se nota como A(x). 

Dada la definition del conjunto de variables que maneja un controlador borroso, sus 
dominios de valores numericos y cualitativos y las funciones de distribucion de posibilidad 
de estos ultimos, una base de n reglas equivalente a la anterior, expresada mediante las 
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funciones de distribution de posibilidad que definen a las etiquetas linguisticas o conjuntos 
borrosos, con la nueva notacion seria la siguiente: 

Ri: Ai(s), Bi(e), Ci(Ae) —> Di(v) 


Ri: Ai(s), Bi(e), Ci(Ae) -»■ D,(v) 


R n : A n (s), B n (e), C n (Ae) —► D n (v) 

Dado un controlador borroso asi dclinido, el problema de inferencia que se plantea es el 
siguiente: en un instante dado, las variables de estado del controlador borroso (s, e y Ae) 
toman los valores s = A’,e= B’ y Ae = C’, cuyas funciones de distribucion de posibilidad 
son, respectivamente, A’(s), B’(e) y C’(Ae); calcular el valor que debe tomar v, notado como 
D’, cuya distribucion de posibilidad es D’(v), para mantener (o devolver) al sistema de 
control a su situacion de equilibrio. 

Este es un problema que se puede resolver mediante la regia composicional de 
inferencia, que, se recuerda, se enuncia como: 

pB-(y) = Sup x {T (pA’(x), p A ->ii(x,y))} 

En el caso de un controlador borroso: 

• La variable y es v. La distribucion de posibilidad a calcular, pB-(y), expresada como 
D’(v), es la accion a realizar para mantener al sistema en equilibro. 

• x es un vector tridimensional: s, e y Ae, donde la distribucion de posibilidad del hecho 
Pa’(x) se calcula mediante la t-norma T (cumple las propiedades asociativa y 
conmutativa) de las funciones de distribucion de posibilidad de los valores que toman las 
variables de estado en un momento dado: A’(s), B’(e) y C’ (Ae). El resultado de esta t- 
nonna representa el resultado de la conjuncion entre A’,B’y C’(A’ A B’ A C’), que son los 
valores cualitativos que toman las variables de estado en un momento dado (el hecho que 
se produce o estado actual del sistema). 

• pA^B(x,y) es Ri: Ai(s) A Bi(e) A Ci(Ae) —> Di(v), para la regia i-esima de controlador 
borroso. 

Por tanto, la regia composicional de inferencia aplicada a la regia i-esima del controlador 
borroso queda enunciada como: 

Di’ (v) = Sup s ,e,Ae {T ((A’(s) A B’(e) A C’ (Ae)),(Ai(s) A Bi(e) A Ci(Ae) D;(v)))} 

Los controladores borrosos emplean la funcion t-norma del minimo: T(x, y)= min (x, y), V 
x G3T, con x, y funciones de distribucion de posibilidad; y la funcion de implicacion de 
Mamdani, delinida como: J (x, y) = min (x, y), V x G.T. Estas restricciones permiten 
simplificar la regia composicional de inferencia definida para controladores difusos, ya que 
tanto las conjunciones como las implicaciones se resuelven mediante la funcion min (x, y). 
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Por tanto, sustituyendo las conjunciones e implicacion de la regia composicional de 
inferencia por la funcion min (x, y), queda la siguiente expresion: 

Di’ (v) = Sup s , e ,Ae {min ((min (min (A’(s), B’(e)), C’ (Ae))), min ((min (min (A;(s), B;(e)), 

Ci(Ae)), Di(v)))} 

Aplicando la propiedad asociativa, se reordenan los terminos y queda la siguiente 
expresion: 

D;’ (v) = min (min (min (Sup s {min (A’(s), A;(s))}, Sup e {min (B’(e), Bi(e))}), SupAe {min 

(C’(Ae), Ci(Ae))}), D,(v)). 

La expresion Sup s {min (A’(s), A;(s))} calcula el valor supremo de posibilidad de la 
funcion de distribucion de posibilidad que representa a la conjuncion entre la etiqueta 
lingiiistica Ai, correspondiente a primera clausula del antecedente de la regia i-esima; y el 
hecho, situacion que ocurre en un momenta dado o estado actual del sistema para esa misma 
variable s, A’. La figura 4.2 muestra graficamcntc como calcularlo y el resultado, que es el 
mayor valor de posibilidad (un numero) que alcanza la funcion de distribucion de posibilidad 
que representa la conjuncion entre A’ y Ai, cuyas funciones de distribucion de posibilidad, 
A’(s) y Ai(s), estan expresadas en funcion de la misma variable s. Tal y como muestra la 
figura 4.2, este valor de posibilidad supremo se nota como NAi,;: nivel de ajuste de la primera 
clausula del antecedente de la regia i-esima A;(s) con la funcion de distribucion de posibilidad 
del hecho A’(s). Si la interseccion entre A’ y Ai es vacia, entonces Sup s {min (A’(s), Ai(s))} = 
0, por lo que NAi,i = 0. Esto significa que la clausula Ai(s) no ajusta con A’(s), por lo que la 
regia R; no se dispara. Por el contrario, si Sup s {min (A’(s), A;(s))} = 1, entonces nos 
encontramos en la situacion de un ajuste nitido. Los valores intermedios de NAi,i indican 
diferentes gradaciones en el ajuste de la clausula del antecedente de la regia con la 
distribucion de posibilidad de su correspondiente hecho, que es una caracteristica propia de la 
logica borrosa. 



Figura 4.2. Valor supremo de posibilidad de la funcion de distribucion que representa a 

la conjuncion entre Ai y A’. 

De la misma manera que Sup s {min (A’(s), A;(s))}= NAi.i, se tiene que Sup e {min 
(B’(e), Bi(e))} = NA23 y SupAe {min (C’(Ae), Ci(Ae))} = NA34, y asi sucesivamente con 
cuantas clausulas existan en los antecedentes de las reglas borrosas de la base de 
conocimiento del controlador. Por tanto: 
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Di’ (v) = min (min (min (NAi,i, NA 2 ,i), NA3,i), Di(v)). 

min (min (NAi,i, NA2.1), NA34) = NA,, que es el nivel de ajuste total del antecedente de la 
regia i-esima con el estado actual del sistema, 0 nivel de ajuste de la regia i-esima. Resulta 
logico pensar que si los valores cualitativos de las clausulas de los antecedentes de las reglas 
de la base de conocimiento del controlador borroso estan unidas mediante conjunciones, el 
nivel de ajuste total del antecedente de la regia sea el minimo de los niveles de ajuste de sus 
correspondientes clausulas. Por esa razon, si alguna clausula j del antecedente de una regia i 
no ajusta con su correspondiente hecho (el que depende de la misma variable) y, por tanto, su 
nivel de ajuste Nji = 0, la regia R, no se dispara: su nivel de ajuste es cero: NAi = 0. 

• La expresion para el calculo de la funcion de distribucion de posibilidad del 
consecuente segun la regia i-esima, Di’ (v), que se corresponderia con pBi’(y) en la regia 
composicional de inferencia, queda como: 

D;’ (v) = min (NAi, Di(v)) 

Donde NA; se debe interpretar como una funcion de distribucion de posibilidad que 
devuelve el valor NAi para todo el dominio; min es, en realidad, la funcion t-norma del 
minimo y Di(v) es el consecuente de la regia i-esima (distribucion de posibilidad). La figura 
4.3 muestra graficamcntc este calculo. 



V 


Figura 4.3. Calculo del consecuente Di’ (v) segun la regia i-esima a partir del calculo de 

su nivel de ajuste Ni. 

• Siguiendo el procedimiento anterior de inferencia en controladores borrosos, para cada 
una de las reglas Ri de la base de conocimiento del controlador borroso, se obtiene una 
funcion de distribucion de posibilidad del consecuente que depende de la variable v: 
{Di’ (v), ..., Di’(v), ..., Dn’(v)}. Por tanto, es necesario calcular la funcion de distribucion 
de posibilidad final D’(v), correspondiente a pB’(y) en la regia composicional de 
inferencia. Para ello, se calcula la t-cononna S (cumple las propiedades asociativa y 
conmutativa) de las funciones de distribucion de posibilidad Di’ (v) V i. El resultado de 
esta t-conorma representa el resultado de la disyuncion (union) entre Df, ..., Di’, ..., D n ’, 
es decir: D’ = Ui (D’i), V i; que son los valores cualitativos, desconocidos (no tienen una 
etiqueta linguistica defmida) por no estar declarados en el controlador borroso, 
correspondientes a las funciones de distribucion de posibilidad obtenidas como resultado 
de la aplicacion de los pasos anteriores. 

Los controladores borrosos emplean la funcion t-conorma del maximo: S(x, y)= max (x, 
y)» v x G3T, con x, y funciones de distribucion de posibilidad expresadas en funcion de la 
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misma variable x, que es dual a T(x, y)= min (x, y) con respecto a la negacion N(x) = 1-x, 
con 1 = pi(x) = 1, V xG3T. 


La figura 4.4 muestra el calculo de la funcion de distribucion de posibilidad D’(v) a partir 
de dos distribuciones de posibilidad D,’(v) y Dj’(v), obtenidas como consecuencia del 
proceso de inferencia en controladores borrosos, suponiendo que unicamente se disparan las 
reglas Ri y Rj. Esto es, que el nivel de ajuste de las reglas de la base de conocimiento del 
controlador es cero, excepto para Ri y Rj, cuyos valores son NA, y NAj, respectivamente. 



Una vez calculada la funcion de distribucion de posibilidad D’(v), resultado del proceso de 
in ferencia con controladores borrosos, es necesario llevar a cabo la interpretacion, conversion 
o desborrocificacion de esta distribucion de posibilidad. Tanto D’(v) como su valor 
cualitativo asociado D’ son, en general, desconocidos porque no estan definidos en el 
controlador borroso. Para llevar a cabo el proceso de interpretacion de D’ (v) se pueden 
emplear cualquiera de las tres tecnicas vistas anteriormente, segun se desee una conversion 
de D’ (v) a un valor nitido o a una etiqueta linguistica de entre las pertenecientes al dominio 
de v. En el primer caso, se calcula el centro de gravedad de D’ (v), gD’, y en el segundo, se 
puede optar por el metodo de la distancia o la X cuenta. 

Generalmente, cuando el receptor de la salida proporcionada por un controlador borroso es 
una persona, como es el caso de los sistemas de ayuda a la decision, la presentacion de 
resultados, alarmas o informacion en pantalla, suele ser conveniente realizar el proceso de 
interpretacion de D’ (v) como una etiqueta linguistica. Las personas somos muy habiles en el 
manejo del conocimiento con imprecision y nos proporciona una informacion mas util el uso 
de etiquetas linguisticas en lugar de valores nitidos precisos. Por ejemplo, en el caso de un 
sistema de calefaccion mediante radiadores, resulta mas util para una persona que el 
controlador borroso informe, mediante la pantalla del termostato, que el agua que circula por 
los radiadores esta “muy caliente” (etiqueta linguistica definida mediante su correspondiente 
funcion de distribucion de posibilidad) o, incluso muestre una serial de alerta, para saber que 
no hay que tocar los radiadores por riesgo de quemaduras. Esta etiqueta linguistica, obtenida 
por el metodo de la distancia o £ C uenta, proporciona mas informacion (a pesar de ser 
imprecisa) que saber el valor nitido exacto; por ejemplo, que la temperatura del agua esta a 
63,5 °C (a partir del calculo del centro de gravedad). En este ultimo caso, podriamos no saber, 
o dudar, acerca de si los radiadores estan calientes. 

Por el contrario, cuando el controlador borroso se comunica con un dispositivo electronico, 
generalmente resulta mas conveniente proporcionar valores nitidos. Continuando con el 
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ejemplo del sistema de calefaccion, resulta del todo inapropiado que el controlador borroso 
indique a la caldera que el agua que circula por los radiadores debe estar “muy caliente” 
porque, seguramente, la caldera no sabra interpretar este resultado. Por el contrario, indicar a 
la carrera que el agua que circula por los radiadores debe estar a 63,5 °C es una accion que 
deberla poder implementar. 


a. Procedimiento general de inferencia en controladores 

Para mantener a un sistema en equilibrio, se fonnula un modelo de control borroso o 
difuso de forma que en cada momenta, se determine el valor que debe tomar una variable v. 
Por tanto, el controlador difuso debe actuar sobre la variable v, la cual viene determinada en 
funcion de las variables de estado s, e y Ae. Para cada una de estas variables, se define su 
significado, dominio de valores numerico y conjunto de valores cualitativos. En este ultimo 
caso, se tiene que s = {Ai, ..., Ai, ..., A n }, e = {Bi, ..., Bi, ..., B n }, Ae = {Ci, ..., Ci, ..., C n } y 
v = {Di, ..., Di, ..., D n }, con posibles repeticiones en todos los casos, es decir, el cardinal de 
estos conjuntos no tiene porque ser n. As! mismo, para cada valor cualitativo del dominio de 
cada una de las variables, se definen sus correspondientes funciones de distribucion de 
posibilidad: Ai(s), ..., A;(s), ..., A n (s); Bi(e), ..., B;(e), ..., B n (e); Ci(Ae), ..., Ci(Ae), ..., 
Cn(Ae); y Di(v), ..., Di(v), ..., D„(v). 

Se dispone de la siguiente base de conocimiento en cualquiera de sus dos notaciones: 

Ri : Ai(s), Bi(e), Ci(Ae) —> Di(v); Ri : s = Ai, e = Bi, Ae = Ci —* v = Di 

Rii A;(s), Bi(e), Ci(Ae) —> D;(v); Ri: s = Ai, e = Bi, Ae = Ci —* v = D; 

Rn: An(s), B n (e), Cn(Ae) > Dn(v)j Rn: s An, e Bn, Ae Cn * V Dn 

En un momenta dado, se sabe que el estado actual del sistema (hecho) es: s = A’, e = B’ y Ae 
= C\ Calcular el valor que debe tomar v para mantener al sistema en equilibrio. 

Para poder calcular el valor que debe tomar v para mantener al sistema en equilibro es 
necesario aplicar el procedimiento de inferencia en controladores difusos, el cual es un caso 
particular de la regia composicional de inferencia en donde T(x, y)= min (x, y); S(x, y)= max 
(x, y), N(x)= 1 - x y J(x, y)= min (x, y); V x G3T, con x, y funciones de distribucion de 
posibilidad expresadas en funcion de la misma variable x y 1 = pi(x) = 1 V x G SC. Asi 

mismo, las clausulas de los antecedentes de las reglas estan unidas mediante conjunciones 
(uniones), que se calculan mediante la t-nonna de las correspondientes funciones de 
distribucion de posibilidad. 
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Este procedimiento general se resumen en los siguientes pasos: 

1. Calculo del nivel de ajuste de cada regia NAi. 

NAi = min {Sup s [T (Aj(s), A’(s))], min (Sup e [T (Bi(e), B’(e))], Sup A e [T (Ci(Ae), C’(Ae))])}, 
siendo T la t-norma T(x, y)= min (x, y). 

Para cada regia Ri, se calcula la t-nonna de las funciones de distribucion de posibilidad 
correspondientes a las etiquetas lingiilsticas de cada una clausulas del antecedente de la regia 
(Ai, B; y Ci, en este caso) con las funciones de distribucion de posibilidad correspondientes a 
las etiquetas lingiilsticas del estado actual del sistema (A’, B’ y C’). Evidentemente, estas 
operaciones t-norma se aplican entre funciones de distribucion de posibilidad dependientes de 
la misma variable: T (Ai(s), Ai’(s)), T (Bi(e), Bi’(e)) y T (Ci(Ae), Ci’(Ae)). Cada una de estas 
operaciones devuelve como resultado otra funcion de distribucion de posibilidad que 
representa la conjuncion o interseccion entre las correspondientes etiquetas linguisticas: Ai A 
A;’, Bi A B;’ y Ci A Ci’. Posterionnente, se calcula el valor supremo de posibilidad de cada una 
de las distribuciones de posibilidad obtenidas como resultado del calculo de la t-norma. 
Finalmente, el nivel de ajuste de la regia Ri es el mlnimo de los supremos calculados 
anteriormente. 

A efectos practicos, en lugar de calcular los niveles de ajuste de todas las reglas de la base de 
conocimiento del controlador borroso, es posible realizar una simplificacion, calculando 
unicamente los niveles de ajuste de aquellas reglas de la base de conocimiento que se 
disparan porque su nivel de ajuste es mayor que cero: NAi > 0. Dado que las etiquetas 
linguisticas de las clausulas de los antecedentes de las reglas estan unidas mediante 
conjunciones, en cuanto se observe, a traves de las correspondientes distribuciones de 
posibilidad, que la interseccion entre cualquier etiqueta lingulstica del antecedente de una 
regia con la etiqueta del hecho correspondiente (aquellas cuyas funciones de distribucion de 
posibilidad dependen de las misma variable) es vacla, la regia no se dispara. Es decir, que si 
se cumple que T (Ai(s), Ai’(s)) = 0 6 T (Bi(e), Bi’(e)) = 0 6 T (Ci(Ae), Ci’(Ae)) = 0, con 0 = 
po(x) = 0 V x G 3C\ entonces la regia Ri no se dispara, siendo NAi = 0. 

2. Calculo de los consecuentes de las reglas que se disparan, NAi> 0. 

Para cada una de las reglas Ri cuyo NAi > 0, se calcula la funcion de distribucion de 
posibilidad del consecuente de la regia Di’(v) como: 

Di’ (v) = min (NAi, Di(v)) 

Donde NA, se debe interpretar como una funcion de distribucion de posibilidad que devuelve 
el valor NAi para todo el dominio; min es, en realidad, la funcion t-norma del mlnimo y Di(v) 
es el consecuente de la regia i-esima (distribucion de posibilidad). 
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3. Calculo de la funcion de distribucion de posibilidad resultado del proceso de inferencia. 

La funcion de distribucion final del proceso de inferencia D’(v) se calcula mediante la t- 
norina (cumple las propiedades asociativa y conmutativa) de las funciones de distribucion de 
posibilidad Di’(v) obtenidas en el paso anterior: 

D’ (v) = T (Di’(v), T(D 2 ’(v), T(Di’(v), ... D n ’(v))...) 

Utilizando la notacion en base a etiquetas linguisticas D’ = El resultado es una funcion de 
distribucion de posibilidad cuya etiqueta linguistica D’ = Ui (D’i), V i. 

4. Interpretacion del resultado final. 

Generalmente, tanto D’ corno su correspondiente funcion de distribucion de posibilidad D’(v) 
son desconocidas por no estar definidas en el controlador borroso. Por tanto, es necesario 
llevar a cabo un proceso de interpretacion o conversion de D’(v) a un numero, mediante el 
calculo de su centro de gravedad, o a una etiqueta linguistica definida en el controlador 
borroso (Di, ..., Di, ..., D n ), mediante los metodos de la distancia o Xcuenta. 

Para el calculo del centro de gravedad de D’(v), se eligen, en primer lugar, puntos de 
discretizacion {vi, v 2 , ..., v m } a intervalos regulares. Posteriormente, se calcula el centro de 
gravedad gD- como: 

m 

S v r D '( v j) 

X d '( Vj ) 

j=l 

Para el calculo de la etiqueta linguistica D; cuya funcion de distribucion de posibilidad se 
parece mas D’(v) mediante el metodo de la distancia, se calculan las distancias de D’(v) a 
cada una de las distribuciones Di (v), ..., Di (v), ..., D n (v) conocidas mediante la siguiente 
formula: 

d(D'(v),D,(v)) = ^/oc-(g D .-g D .) 2 + (3 ■ (A D , - A D . ) 2 

Siendo a y /? valores de ponderacion, a+/3= 1 y a, p [0,1]; gD-, gDi los centros de gravedad y 
Ad’, Adi las areas de D’(v) y Di(v), respectivamente, aquella D, cuya funcion de distribucion 
de posibilidad, Di(v), que minimiza d (D’(v),Di(v)) es la etiqueta linguistica que se devuelve 
como resultado: Di / V j = l..n, d (D’ (v), Di(v)) < d (D’ (v), Dj (v)) 

Por ultimo, para el calculo de la etiqueta linguistica Di cuya funcion de distribucion de 
posibilidad se parece mas D’(v) mediante el metodo de la £ C uenta, se calcula, para cada una 
de las distribuciones Di (v), ..., D;(v), ..., D n (v) conocidas : 

V (r' a a ^ 

X-JV B ')= y fA V 

Aquella Di cuya distribucion de posibilidad D; (v) maximiza Xcuenta (D; / D’) es la etiqueta 
linguistica que se devuelve como resultado: Di / V j = 1 ..n, Y, cuenta (Di/D’) ^ Ycuenta (Dj/D’). 
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A lo largo de todo el desarrollo del proceso de inferencia en controladores borrosos, se 
ha supuesto que las reglas de la base de conocimiento tienen tres clausulas, dependientes de 
las variables s, e y Ae. Evidentemente, pueden tener cualquier numero de clausulas, variando 
unicamente el numero de terminos que intervienen en el calculo de los niveles de ajuste de 
las reglas: NAi = min {Sup s [T (Ai(s), A’(s))], min (Sup e [T (Bi(e), B’(e))], SupAe [T (Ci(Ae), 
C’(Ae))])}. Un caso particular se produce cuando alguna de las reglas (Ri) tiene una unica 
clausula, por ejemplo, dependiente de la variable s. En este caso, el nivel de ajuste es el valor 
supremo de la funcion de distribucion de posibilidad resultante de aplicar la t-norma entre la 
funcion de posibilidad A;(s) y la correspondiente al hecho que depende de la misma variable 
A’(s). Por tanto, NAi = Sup s [T (A;(s), A’(s))]. 

b. Ejemplo de inferencia en controladores 

Para controlar el sistema de calefaccion de una vivienda, se ha disenado un modelo de 
control borroso, de forma que, en cada momento, se detennine cual debe ser la temperatura 
del agua que circula por los radiadores con el fin de mantener el ambiente a la temperatura 
deseada por el usuario. 

El sistema de calefaccion consta de los siguientes componentes: Un tennostato 
electronico. Este es el aparato donde va empotrado el sistema de control borroso para 
controlar la temperatura ambiente. El tennostato posee un reloj y un programador que 
pennite al usuario establecer diferentes temperaturas ambiente a lo largo del dla. Una caldera 
de gas, encargada de calentar el agua que circula por los radiadores bajo mandato del 
tennostato. Los radiadores, distribuidos por la vivienda. Finalmente, los sensores de 
temperatura, conectados al tennostato para indicar la temperatura ambiente actual. 

El controlador difuso debe actuar sobre la variable Temperatura del Agua que circula 
por los radiadores (TA) para controlar la temperatura ambiente. Esta variable viene 
detenninada en funcion de la Diferencia entre la Temperatura ambiente actual y la 
Temperatura definida por el usuario en el tennostato (DT), la variacion de esta diferencia de 
temperaturas (dDT) y la Hora actual (H), tomada del reloj del tennostato. Por tanto, se 
manejan las siguientes variables: 

• TA: Temperatura del Agua que circula por los radiadores, con valores posibles: Fria 
(F), Templada (T), Caliente (C) y Muy Caliente (MC). Su intervalo de valores es 
(0,100)°C. 

• DT: Diferencia de Temperaturas entre la que existe actualmente en el ambiente y la 
definida por el usuario en el tennostato, con posibles valores: Por Debajo (PD), 
Exacto (E) y Por Encima (PE). Su intervalo de valores [-50,50] °C. Los valores 
negativos indican que la temperatura ambiente esta por debajo de la deseada por el 
usuario, mientras que los valores positivos indican que la temperatura ambiente esta 
por encima de la deseada por el usuario. 
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• dDT: Variacion de la Diferencia de Temperaturas, medido en grados centlgrados por 
hora (°C/h.), con valores posibles: Negativa (N), cero (Z) y Positiva (P). Su intervalo 
de valores es [-5,5]°C/h. 

• H: la Hora actual, con posibles valores: Manana (Ma), Tarde (Ta), Noche (No) y 
Madrugada (M). Su intervalo de valores es [0,24) horas. 

Los valores cualitativos que pueden tomar las variables que maneja el controlador 
borroso tienen asociadas las funciones de posibilidad mostradas en la figura 4.5: 





Figura 4.5. Funciones de distribucion de posibilidad asociadas a los valores cualitativos 

de TA, DT, dDT y H. 

Para relacionar las variables, se dispone de la siguiente base de reglas: 


Regia 

DT 

dDT 

H 

* TA 

Rl 

— 

— 

Ma 

F 

R2 

PE 

Z 

No 

F 

R3 

E 

P 

No 

F 

R4 

PE 

N 

Ta 

T 


Se pide aplicar el modelo de razonamiento en controladores borrosos para obtener el 
valor numerico (para enviar a la caldera) y cualitativo (para mostrar en la pantalla del 
tennostato) de TA (temperatura del agua) que debe circular por los radiadores, sabiendo que 
es por la rnanana (Ma), que la diferencia de temperaturas es de +2 °C y que la variacion de la 
diferencia de temperaturas es de -2,5 °C. 

SOLUCION: 


1. Calculo del nivel de ajuste de cada regia NAi. 

En primer lugar es necesario calcular el nivel de ajuste NAi del estado actual del sistema con 
el antecedente de cada una de las reglas Ri que forman la base de conocimiento del 
controlador borroso. En lugar de calcular los niveles de ajuste de todas las reglas de la base 
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de conocimiento del controlador borroso, es posible realizar una simplificacion, consistente 
en calcular unicamente aquellos NAi>0. Para ello, es necesario conocer, a traves de las 
correspondientes distribuciones de posibilidad, que la interseccion entre cualquier etiqueta 
lingulstica del antecedente de una regia con la etiqueta del hecho correspondiente (aquellas 
cuyas funciones de distribucion de posibilidad dependen de las misma variable) no es vacla. 
Por tanto, se calcula a continuacion con que etiquetas lingulsticas (conjuntos borrosos) hay 
interseccion no vacla con aquellas que definen el estado actual del sistema de calefaccion, as! 
como el nivel de ajuste entre sus correspondientes funciones de distribucion de posibilidad. 

Se dice en el enunciado que el estado actual del sistema se describe de la siguiente forma: 
H = Ma; DT = +2 y dDT = -2,5. Se sobreentiende que la funcion de distribucion de 
posibilidad correspondiente a DT = +2 es p +2 (DT = +2) = 1; 0 en el resto. De forma similar 
para dDT = -2,5: p- 2,5 (dDT = -2,5) = 1; 0 en el resto. 

• H = Ma tiene interseccion no vacla con madrugada (Ma), Manana (M) y Tarde (T). Se 
calculan entonces los supremos de las t-normas de las distribuciones de posibilidad 
correspondientes al hecho Ma con Ma, M y T (niveles de ajuste): Supu[T (Ma(H), 
Ma(H))] = 1; Sup H [T (Ma(H), M(H))] = 0,5; Sup H [T (Ma(H), T(H))] = 0,5. La figura 4.6 
muestra estos calculos de forma grafica 



Figura 4.6. Calculo de los niveles de ajuste entre las distribuciones de posibilidad 
correspondientes al hecho H = Ma y a aquellas con cuyas etiquetas linguisticas existe 

una interseccion no vacia: Ma, M y T. 

• DT = +2 tiene interseccion no vacla con Por Encima (PE): Supdt[T(+ 2(DT), PE(DT))]=1. 
La figura 4.7 muestra este calculo de forma grafica. 


A 

1 


-50 


PE (DT) 
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- 2 - 10+1 +2 


+50 


DT 


-> 
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Figura 4.7. Calculo del nivel de ajuste entre +2 (DT) y PE (DT). 


• dDT = -2,5 tiene interseccion no vacia con Negativa (N): SupdDT[T(-2,5(dDT), 
N(dDT))]=l. La figura 4.8 muestra este calculo de forma grafica. 
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-5 -2.5 0 +2.5 +5 °C/h. 

Figura 4.8. Calculo del nivel de ajuste entre -2,5 (dDT) y N (dDT). 

Por tanto, solo se disparan aquellas reglas que tengan en sus antecedentes Ma, M o T para la 
variable H; PE para la variable DT y N para dDT; esto es: 

Rl: DT = dDT = -, H = Ma —» TA = F; 

R4: DT = PE, dDT = N, H = Ta —» TA = T; 

Dado que Rl solo tiene una clausula en el antecedente, H = Ma, el nivel de ajuste de la regia 
es el nivel de ajuste de su unica clausula con el hecho: NAi = NAn; siendo NAn = Suph[T 
(Ma(H), Ma(H))] = 1. Por tanto, NAi = 1. 

En el caso de R4, con tres clausulas en el antecedente, NA 4 = min (NA 14 , NA 24 , NA 34 ), donde 
NA 14 = Su P dt[T(+2(DT), PE(DT))]=1; NA 24 = SupdDi[T(-2,5(dDT), N(dDT))]=l; y 
NA 34 = Sup H [T (Ma(H), T(H))] = 0,5. Por tanto, NA 4 = 0,5. 

2. Calculo de los consecuentes de las reglas que se disparan: Ri y R 4 . 

La formula general del calculo de la funcion de distribucion de posibilidad del consecuente 
de una regia en un controlador borroso es la siguiente: 

D;’ (v) = min (NAi, Di(v)) 

Donde NAi se debe interpretar como una funcion de distribucion de posibilidad que devuelve 
el valor NAi para todo el dominio; min es, en realidad, la funcion t-norma del minimo y Di(v) 
es el consecuente de la regia i-esima (distribucion de posibilidad). 

En el caso de este ejercicio, es necesario calcular DT(TA) y D 4 ’(TA), con Di’(TA) = min 
(NAi, F(TA)); y D 4 ’(TA) = min (NA4, T(TA)). Dado que NAi = 1, interpretado como pi(x) = 
1 V x G 3C, es el elemento neutro de la t-norma, en este caso la t-norma del minimo, se tiene 

que Di’(TA) = F(TA), la cual se muestra en la figura 4.5 en color azul. La figura 4.9 muestra 
graficamente el calculo de la funcion de distribucion de posibilidad D 4 5 (TA). 



Figura 4.9. Calculo de la funcion de distribucion de posibilidad Df (TA). 
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3. Calculo de la funcion de distribucion de posibilidad resultado del proceso de inferencia. 

La funcion de distribucion final, D’ (TA), se calcula mediante la t-conorma del maximo de las 
funciones de posibilidad Di’(TA) y D 4 ’(TA): D’(TA) = S (Di’(TA) y D/t’(TA)). La figura 4.10 
muestra graficamente el calculo de D’(TA). 

F (TA) 



Figura 4.10. Calculo de la funcion de distribucion de posibilidad D’ (TA). 


4. Interpretacion del resultado final: D’ (TA). 

Dado que la funcion de distribucion de posibilidad D’(TA) no esta definida en el sistema de 
control borroso propuesto en el enunciado de este ejercicio, es necesario llevar a cabo un 
proceso de interpretacion o conversion de D’(TA). En este caso, se solicita tanto el valor 
numerico (para enviarlo a la caldera del sistema de calefaccion), como la etiqueta linguistica, 
para mostrarla en la pantalla del termostato. 

El valor numerico de TA para enviarlo a la caldera y que, de esa forma, el sistema a controlar 
vuelva a su estado de equilibrio, se lleva a cabo mediante el calculo del centro de gravedad de 
D’(TA), tomando los puntos de discretizacion LA a intervalos de 10°C: 


§D' = 


n 

^TAj DfTA;) 

i=l_ 

i>( TA i) 


10-1 + 20-1 +30-0,5+ 40-0,5+ 50-0,5 

1 + 1 + 0,5 + 0,5 + 0,5 


25,71 


Por tanto, es necesario enviar a la caldera una orden de calentar el agua que circula por los 
radiadores a 25,71°C para devolver el sistema de calefaccion a su estado de equilibrio. 


El valor cualitativo o etiqueta linguistica de las definidas en el controlador borroso cuya 
funcion de distribucion de posibilidad se parece mas a D’(TA), con el fin de enviarlo a la 
pantalla del termostato como informacion al usuario, se calcula mediante la £cuenta. 

Las etiquetas linguisticas conocidas, con las que hay que comparar D’, son F, T, C y MC, 
cuyas funciones de distribucion de posibilidad se muestran en la figura 4.5. En primer lugar 
es necesario calcular la £cuenta de estas funciones de distribucion de posibilidad, para lo que 
se toman puntos de discretizacion a intervalos de 10 °C: 

• Icuenta(F)= 1 + 1 + 0,5 = 2,5. ^cuenta (T) = 0,5 + 1 + 0,5 = 2. 

* y^.cuenta (C) = 0,5 + 1 + 0,5 = 2. Xcuenta (MC) = 0,5 + 1 + 1 = 2,5. 
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Posteriormente, utilizando los mismos puntos de discretizacion, se calculan las Xcuenta de las 
funciones de distribucion de posibilidad que representan la interseccion entre cada una de las 
etiquetas lingulsticas conocidas (F, T, C, MC) con la obtenida tras el proceso de inferencia en 
controladores borrosos, D’. La llgura 4.11 muestra graficamente las funciones de distribucion 
de posibilidad que representan el resultado de estas intersecciones. 



D' (TA) 


A 


S(D'(TA), C (TA)) 


TA 

H> 


Figura 4.11. Funciones de distribucion de posibilidad que representan el resultado de la 

interseccion de D’ con F, T, C y MC. 

* Xcuenta (D’ A F) = 1 + 1 + 0,5 = 2,5. 

* Xcuenta (D’ A T) = 0,5 + 0,5 + 0,5 = 1,5. 


* Xcuenta(D’ A C) = 0,5. 

• Xcuenta (D’ A MC) = 0, dado que la interseccion entre D’ Y MC es vacla. 


Por ultimo, se calcula la £ cuenta (F / D ), ^ cuenta (T/D )? y.cuenta (C/D ) y y.cuenta (MC / D’) 
como: 


X (Aj/D') = 

cuenta 1 


y (d ,a a,) 

cuenta 1 


A. 


Donde Ai representa a F, T, C y MC, las etiquetas lingulsticas cuyas funciones de distribucion 
de posibilidad conocidas, definidas en la figura 4.5. 


* Xcuenta(F/D’) = 2,5/2,5= 1 

* Icuenta(T/D’)= 1,5/2 

* Icuenta (C / D’) = 0,5/2 

* y.cuenta (MC / D’) = 0 


La etiqueta linguistica Ai que maximiza ^cuenta (Ai / D’) es F, por tanto la pantalla del 
tennostato del sistema de calefaccion debe indicar que el agua que circula por los radiadores 
esta frla. 
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